Coupled-channel dynamics in the reactions ttN — > ttN, rjN, KA, KH 
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Elastic nN scattering and the world data of the family of reactions tt p 



rjn, 



K°A, K°Z°, 



K + Tj~ , and ir + p — > K + Y, + are described simultaneously in an analytic, unitary, coupled-channel 
approach. SU(3) flavor symmetry is used to relate both the t- and the u- channel exchanges that drive 
the meson-baryon interaction in the different channels. Angular distributions, polarizations, and 
spin-rotation parameters are compared with available experimental data. Partial-wave amplitudes 
are determined and the resonance content is extracted from the analytic continuation, including 
resonance positions and branching ratios, and possible sources of uncertainties are discussed. The 
results provide the final-state interactions for the ongoing analysis of photo- and electroproduction 
data. 

PACS numbers: 11.30.Hv, 11.80.Gw, 13.75. Gx, 14.20.Gk, 24.10.Eq 



I. INTRODUCTION 



To gain insight into the non-pcrturbativc sector of 
Quantum Chromodynamics (QCD) the knowledge of the 
excited hadron spectrum is essential, providing the con- 
nection between experiment and QCD. An ab-initio ap- 
proach to resonance physics is provided by lattice QCD 
gauge simulations. Spectra of excited baryons observed 
on the lattice [il-flTI] hold the promise to explain the rich 
phenomenology in the first, second, and third resonance 
regions. Once quark masses drop towards the physical 
limit and finite volume effects, increasingly important for 
lighter quark masses, are under full control, comparison 
to experimental data will be possible. Dynamical Dyson- 
Schwingcr approaches [PZHTq and quark models |16H25j 
also rely on a comparison with the extracted resonance 
spectrum. A complementary picture for some resonances 
is p rovided by the unitarization of chiral interactions [2(1 - 
l45j . using directly the hadronic degrees of freedom and, 
as such approaches are formulated for at most a few par- 
tial waves, one again needs a data analysis in terms of 
partial waves to compare to. 

The reliable extraction of amplitudes and of the res- 
onance spectrum from data, therefore, is a prerequi- 
site for comparing theoretical approaches to experiment. 
Most resonances have been identified through elastic irN 
scattering in the past [4^l- l48f and to this day by the 
GWU/SA1D analysis group |49m6l |. A new experimen- 
tal window has opened through the recent high-precision 
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photon-beam facilities, e.g., at ELS A, GRAAL, JLab, 
SPring-8 and MAMI, allowing for resonance identifica- 
tion and the determination of their electromagnetic prop- 
erties mm. 

Combining different reactions for resonance extrac- 
tion enables one to determine those states which couple 
only weakly to TtN. The simultaneous analysis of dif- 
ferent final states of pion- and photon-induced reactions 
in coupled-channel approaches is, thus, the method of 
choice to study reaction dynamics. Dynamical coupled- 
channel (DCC) models |59l - [86| provide a sophisticated 
tool for analyzing excited baryons and extracting reso- 
nance parameters as they obey theoretical constraints of 
the iS-matrix such as analyticity. 

The DCC model developed and employed in this study 
(Jiilich2012 model) is based on an approach pursued over 
many years [59l - l66l | . The scattering amplitude is obtained 
as the solution of a Lippmann-Schwinger-type scatter- 
ing equation defined by time-ordered perturbation the- 
ory (TOPT) which automatically implements unitarity 
constraints. Also, the important ttttN channels are in- 
cluded dynamically in the sense that the 7T7r and ttN 
subsystems match the respective phase shifts. Based on 
effective Lagrangians, the non-resonant interactions are 
given by t- and u-channel exchanges of known mesons 
and baryons, while s-channel processes can be consid- 
ered as bare resonances. The explicit treatment of the 
background in terms of t- and w-channcl diagrams in- 
troduces strong correlations between the different partial 
waves and generates a non-trivial energy and angular de- 
pendence of the observables. Analyticity is respected by 
including dispersive contributions of intermediate states, 
as well as the correct structure of branch points [8714891 ] . 
Thus, a reliable determination of resonance properties 
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given in terms of pole positions and residues of the scat- 
tering amplitude in the complex energy plane is possible. 

The extension of the Jiilich model to photoproduction 
within a gauge-invariant approach [67H69I ] has been ac- 
complished recently [70l |. There are a number of stud- 
ies of photon-induced reactions within phenomenological 
dynamical models 0, Izl HI HI The vast 

majority of them, however, do not satisfy the gauge in- 
variance condition as dictated by the generalized Ward- 
Takahashi identity 1671 [9911 . E xceptions to this are, e.g., 
the works of Refs. [Mlslllol. 

The adaptation of DCC models to the finite volume, 
to allow for the prediction of lattice levels and the cal- 
culation of finite volume corrections, was pioneered in 
Ref. [l OOj ] - In principle, such extensions of hadronic ap- 
proaches allow for the anal ysis of ex perimental and lattice 
data on the same footing [Toil - fTo3| . 

In Ref. (6(| , the Jiilich model is extended to the kaon- 
hyperon sector. SU(3) flavor symmetry is exploited to 
relate different coupling constants; SU(3) symmetry is 
broken e.g. by physical masses but also by different cut- 
offs in the form factors of the vertices. The latter absorb 
other sources of SU(3) breaking that are not systemat- 
ically addressed as possible in chiral perturbation the- 
ory. A resonance analysis of the isospin 7 = 3/2 sector 
was performed by fitting simultaneously the elastic ttN 
partial- wave amplitudes of the GWU/SAID analysis 
and ir + p — > K + T, + differential cross se ctions, polariza- 
tions, and spin- rotation data [I04| - |l06| . To achieve a 
good description of the A + £ + data, the inclusion of res- 
onances with total spin up to J = 7/2 turned out to be 
necessary. 

In the present study, we extend this resonance analysis 
to the isospin 7=1/2 sector, considering the world data 
of the set of reactions n~p ijn, K°A, A°£°, K + E~ , 
and ir + p — > K + Y, + , together with irN — > irN scattering. 
Besides the ttN, rjN, KA and A£ channels, the approach 
includes three effective ttttN channels, namely ttA, aN 
and pN. The considered energy range has been extended 
beyond 2 GeV and resonances up to J — 9/2 are included 
now. For older coupled-channel a nalyses of pion-induced 
inelastic reactions see, e.g., Refs. [lOMQl. 

The present effort is a first step towards a global 
analysis of pion- and photon-induced production of nN, 
ijN, KA and A£. For example, photoproduction data 
with unprecedented accuracy were measured recently at 
ELSA, MAMI, JLab and other facilities for the KY 
final states 110l - ll27| . the latter being also acce s sible 
through nu cleon-nuc leon collisions at COSY [l28l Il29j | 
or HADES ptlJl^ . New data for pion-induced reac- 
tions are expected by experi ments such as EPECUR at 
ITEP for the KA final state fl3l a nd E19 at J-PARC 
for the charged KY, final states [l34j . 

This paper is organized as follows: An overview of the 
formalism is given in Sec.|Hl including a new renormaliza- 
tion scheme in Sec. Ill Bl Details on the data base and the 
fitting techniques are documented in Sec. IIII1 together 
with the comparison of the present results to data. In 




FIG. 1: (Off-sh ell) pion exchange as a consequence of three- 
body unitarity [1551 ] . 



particular, results of two fit scenarios (A and B) are pre- 
sented which differ predominantly in the strength of the 
coupling of the ttN system to the effective irirN channels, 
notably to aN. The description of observables achieved 
in the various reaction channels is discussed in detail. 
Extracted resonance pole properties are presented and 
discussed in Sec. IIVI The technical details and values of 
fit parameters shown in the Appendices ensure the repro- 
ducibility of the results. 



II. FORMALISM 

Theoretical requirements on the scattering amplitude 
include two-body unitarity which is realized in most of 
the current approaches to dat a analysis , including those 
formulated with the A'-matrix p5Wl49j . Dispersive con- 
tributions of the two-body intermediate states are ne- 
glected in A-matrix approaches. However, there are cal- 
culations in which dispersive parts are in cluded in a fac- 
torized form, e.g. CMB-type approaches jl5Cll4l54l | or the 
Chew-Mandelstam formalism used in the GWU/SAID 
analysis [49l - l56j . 

In addition to two-body unitarity, three-body unitarity 
is an important issue when analyzing the baryon spec- 
trum due to the known large inelasticities caused by the 
ttttN channel. In the current approach, the irirN interme- 
diate states are parameterized as ttA, aN, and pN chan- 
nels. The properties of the unstable particles in these 
channels are matched to the corresponding phase shifts, 
i.e. the A fits the P33 partial wave in ttN scattering, 
and the a and p mesons fit the 7T7t phase shifts for the 
corresponding quantum numbers J62|, [65| . It should be 
stressed that the current approach provides the correct 
analytic structure, including branch points in the com- 
plex plane from three-particle intermediate states. In 
Ref. [89| it was shown that the latter are indeed of prac- 
tical relevance, to avoid erroneous resonance signals. 

For three-body unitarity, the consistent inclusion of 
three-body cuts is essential [155] as illustrated in Fig. [T] 
In particular, three-body unitarity dictates the inclusion 
of the discontinuities induced by ^-channel exchanges si- 
multaneously with the discontinuities arising from three- 
body states propagating in the s-channel. Using disper- 
sion relations, the form of the t-channel and s-channel 
propagators can be calculated from these discontinu- 
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ities [15a ]. It should be noted that this implies the full 
off-shell treatment of the exchange potential. In addition 
to t-channel exchanges, it is always possible to add more 
interactions as long as they do not induce cuts, such as 
contact terms. 

Thus, meson exchange arises naturally from require- 
ments of the S'-matrix. While pion exchange is included 
in the transitions potentials (cf. Figs. [5] and [31) it should 
be noted that full three-body unitarity requires addi- 
tional diagrams that are beyond the scope of this work. 
Those comprise, e.g., the dressing of the nucleon in the 
■nN s-channel propagator, as induced by the nucleon ex- 
change term in 7riV — » t tN . This contribution is known 
to be small, though [155} . Baryon-exchange diagrams are 
also mandatory for the approximation of the left-hand 
cut on the amplitude level. The extension to channels 
involving particles with strangeness then requires addi- 
tional diagrams. 

All resulting t- and w-channel processes of the current 
approach are shown in Figs. EH and [31 For the transitions 
MB -> MB, with the J p = 0~ octet mesons M and the 
J p = l/2 + octet baryons B, we have included the full set 
of diagrams required by SU(3) symmetry. This applies to 
the classes of diagrams of the exchange of J p = 1 ~~ vector 
mesons, J p = 1/2+ baryons, and J p = 3/2+ baryons. In 
particular, non-diagonal transitions such as r/N — > KY 
are considered. 

At this point, we would like to stress that in the current 
approach the ttN interactions with the cr(600) and p(770) 
quantum numbers in the f-channel, as shown in Fig. 21 
are determined from analytically continued NN — > tttt 
data, using dispersive techniques and crossing symme- 
try (see Ref. j59[ and references therein). Note that 
for low energies, where the ititN channel plays no sig- 
nificant role, more systematic schemes using dispersive 
techniques for th e amplitu de extraction have been de- 
veloped recently [l56l Il57l | . Chiral perturbation theory 
also allows for a model- i ndep enden t understanding of the 
near-threshold behavior [I58l4l62| . 

The interaction potentials V, given by t- and u-channel 
exchanges of hadrons as shown in Figs. [H and EH ar e de- 
rived from the Wcss-Zumino interaction [IM m$ plus 
additional terms. These Lagrangians, in combination 
with the form factors, determine also the off-shell be- 
havior of the amplitude. For analytic expressions of the 
exchange diagrams, see Appendix [B] The Lagrangians 
and vertices for the s-channel exchanges are listed in Ap- 
pendix B of Ref. (6|| ; see also Appendix O of this study. 

In the present approach, an integral equation is solved 
to obtain the T-matrix based on the exchange processes. 
From the resulting T-matrix, the observables in pion- 
induced reactions can be calculated [66| (for photon- 
induced reactions, see Refs. [68l47fjl ]l. The scattering 
equation in the center-of-mass (cm.) frame, for a given 



partial wave, reads [62M6 1 
T flI/ (p",p',z) = V li „(p",p',z) 



V^(p",P,z)T KV (p,p',z) 
z - E a (p) - E b (p) + ie 



(1) 



where p" = \p"\ (p' = \p'\) is the modulus of the outgo- 
ing (incoming) three-momentum that may be on- or off- 
shell, z = i/s = E is the scattering energy, and p, v, n 
are channel indices. Note also that the different possibil- 
ities of the 7rA and pN states to couple to a given J p are 
organized in ad ditional c hannels as li sted in T able IIX1 In 
Eq. ([1]), E a = y/ m 2 + p 2 and Eb = \Jm 2 + p 2 are the en- 
ergies of the intermediate particles a and b in channel n. 
Eq. (TTJ) is formulated in the partial-wave basis, i.e. the 
amplitude only depends on the modulus of the incom- 
ing, outgoing, and intermediate particle momenta. This 
implies a partial-wave decomposition of the exchange po- 
tentials [62|,[63[ (see also Table [IX] in Appendix [B]) . The 
angular dependence of the full T-matrix is provided by 
the Wigner d{\, (^)-functions in the partial-wave decom- 
position [62l . |63| , where $ is the scattering angle and 
A (A') is the helicity of the incoming (outgoing) baryon. 
The denominator in Eq. (TTJ) corresponds to the channels 
with stable particles, irN, r)N, KA, and KY; for the ef- 
fective ttttN channels (ttA, aN, pN), the propagator is 
more complex (62l . [65j . While in the present approach the 
resonances acquire their widths from the solution of the 
scattering equation ([1]), there are also other approaches 
to unst able particles (see, e.g., the complex-mass scheme 
of Ref. [TH; see also R ef. jl66| ). 

It should be stressed that throughout this study we 
use isospin-averaged masses. For the different KY, fi- 
nal states considered, in principle, the threshold energies 
are different which could have quantitative effects close 
to threshold; however, all KY, data available lie consider- 
ably above the KY thresholds and effects from differences 
in those thresholds are negligible. Similar considerations 
hold for electromagnetic effects. For elastic ttN scatter- 
ing, we do not fit to data anyway but to the GWU/SAID 
amplitudes [54j (cf. Sec. IIII C[) which include corrections 
for Coulomb effects. 

Having outlined the approach, it is obvious that hadron 
exchange connects all partial waves of the scattering am- 
plitude, putting strong constraints on the form of the 
amplitudes and providing at the same time a background 
from t- and w-channel processes. Together with a min- 
imal set of s-channel processes ("genuine" resonances), 
this allows for a reliable extraction of poles and residues 
from the analytic continuation of the amplitude as de- 
rived in Ref. 1651. 



A. Decomposition of the amplitude 

The scattering amplitude can be decomposed into 
a non-pole and a pole part as derived in this section. In 
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FIG. 2: Hadron exchanges in the t- and it-channel employed in the current approach. Not shown are s-channel processes 
(resonances) . 



the present study, this decomposition merely serves to 
speed up to the numerical evaluation; no physical mean- 
ing is attached to the components of the decomposition. 

The sum of the u- and t-channel diagrams is labeled 
as V NP in the following. Together with the (bare) s- 
channel exchanges V p , they constitute the interaction V 
in Eq. Q, 

V, v = C P + V* = C P + f; ^% (2) 

i— o 1 

with n being the number of bare s-channel states in a 
given partial wave. To simplify the notation, the explicit 
dependence on momenta and energy arc omitted here 



and in the following. In Eq. @, j c ["f a ] are the bare 
resonance creation [annihilation] vertices, indicated by 
the subscript c [a] , of the s-channel states of bare mass 
mf-. The bare vertices for J < 3/2 are derived from 
Lagrangians and given in Appendix B of Ref. (66j . The 
vertex functions for J > 5/2 are given in Eq. (B.3) of 
Ref. [Hi and Eq. ([AT]) in this study. 

The u-channcl diagrams from nuclcon, A(1232), A, S, 
£*(1385), S and S*(1530) exchanges are included with 
physically known coupling strengths, while u-channcl di- 
agrams from other baryonic resonances are neglected. 
Those would introduce additional parameters which are 
difficult to adjust, for the diagrams do not introduce 
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FIG. 3: Continuation of Fig. H 



strong energy depende ncies (for a discussion of u-channcl 
contributions, see Ref. [167I ]). 

The unitarization of 1/ NP leads to the definition of 

the so-called non-pole part of the full T matrix (projected 
to a partial wave). The G K are the propagators of two- 
or three-particle intermediate states [i.e., for the two- 
particle case, G K = (z — E a — + ie)^ 1 according to 
Eq. (HJ]. Note that in Eq. ©, and also Eq. (g]) be- 
low, there is an integration over the momentum of the 
intermediate state, cf. Eq. ([!}, not written explicitly for 



simplicity. 

Contributions from the s-channcl exchanges, the pole 
part T p , can be evaluated from the non-pole part T NP 
given in Eq. For this, we define the quantities 

s u = E^^r^. ( 4 ) 

where T c [r a ] are the dressed resonance creation [anni- 
hilation] vertices and £ is the self-energy. The indices 



6 



t (NN-» mi) 




S (;iN-> ttN) It 



FIG. 4: Exchanges with cr(600) and p(770) quantum numbers 
in the i-channel, determined by using dispersion techniques 
and crossing symmetry 



i, j label the s-channel state in the case of multiple reso- 
nances. For the two-resonance case, the pole part reads 
explicitly [36| 



D = 



T^D- 1 T c l/ , where 

/pa pn \ pc 

Z — — En 

-E 2 i z 



pc 

-E 



12 



(5) 



from which the one-resonance case follows immediately. 

It is easy to show that the full scattering T matrix of 
Eq. ([T]) is given by the sum of pole and non-pole part, 



T — T p 



■T, 



NP 



(6) 



This deco mposition is widely used in the literature, see 
e.g. Refs. 0,[IMG]iI- The non-pole part T NP is some- 
times referred to as background, although the unitariza- 
tion of Eq. may lead to dynamically generated poles in 
T NP . In turn, T p contains constant and other terms be- 
yond the pole contribution of the Laurent expansion [64| ■ 
Consequently, the conclusion was drawn that the clean- 
est separation into a background and a resonance part is 
given by the separation into a singularity-free part and 
the part a_i/(z — z ) that is the leading term in the Lau- 
rent expansion of the scattering amplitude in the complex 
energy plane about the pole position z = zq with residue 
a_i. 

However, in the present study, we use the decompo- 
sition of Eq. ([5]), because the calculation of T p is nu- 
merically much faster than that of T NP . In a fit of 
only s-channel parameters, it is thus convenient to cal- 
culate T NP once and then fit the resonance parameters, 
which only requires the multiple re-evaluation of T p (cf. 
Sec. HUB]). 



B. Renormalization of the nucleon mass and 
coupling 

The nucleon is included in our formalism as an explicit 
s-channel state which is dressed via Eq. ((4]) to acquire 
the correct physical nucleon mass and ttNN coupling. 



This means, to ensure a nucleon pole at the position of 
the physical nucleon mass tun = 93 9 M eV and with the 
physical coupling of J^nn = 0.989 |l70j . one has to im- 
plement a renormalization procedure in the P\\ partial 
wave. 

The renormalization in case of one bare s-channel 
state, coupling to one channel (tvN), is w e ll un derstood 
and has been applied, e.g. in Refs. }l69l Il7l| . For il- 
lustration, we repeat the essential steps in the current 
notation. 

Following Eq. ([5]), the pole part T p in the one- 
resonance case, with an s-channel nucleon exchange of 
bare mass m b N , can be written as 



T p = 



1 1 1 1 



En 



(7) 



where here and in the following the omission of the chan- 
nel index means that we refer to the wN channel. We 
require that T p has a pole at the physical nucleon mass 
m N : 

m b N + En (zq) with z 



(8) 



To derive an expression for m b N that respects the renor- 



malization condition of Eq. 
self-energy about z = zq 

En(z) = Eii(z ) + (z - Zq) 



we expand the nucleon 
dSu(z) 



■o(z- Zo y 



(9) 



and define the reduced quantities En and f°' c taking 
advantage of the fact that the bare coupling ffj can be 
factorized off E and T, 



En 



( rb \ 2 v-ii'cd 
UN) ^11 

-ircd:a,c 



2 r 2 yrcd 
x J-kNN ^11 

^rcdia.c 



2 

X En 



1 = Jni =x}^nn^ 1 =xT 1 ' , (10) 



with x £ 



and the bare irNN coupling constant 



N 



NN 



(11) 



In Eq. (flQ|) . En is the nucleon self-energy calculated with 
the physical nucleon coupling f^NN, instead of the bare 
coupling. The same applies to r"' c . T p from Eq. ([7]) 
reads now 



T p = 



1 



x^lfl 



Z- Zq 



1 



c 2 9,En 



+ O(z-z )° (12) 



with d z := 10 determine tne Dare coupling /j^, x 
is calculated: the physical residue of the nucleon pole in 
the variable z is given by 



JV-s-ttJV 



7i c 7? 



(13) 



where 7"'° are the bare nucleon vertices calculated at 
z = zq with the physical nucleon coupling f^NN instead 
of the bare coupling (cf. Appendix B.l. of Ref. |6o|). 



7! 



. / 3 /ttjvat k 
i \ — 



E 



N 



m N 



8 Trm^ yjE N bj v {E N + m N ) 



(14) 
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where k is the particle momentum in the center-of-mass. 
The physical residue of Eq. ([13")) has to agree with the 
residue of T p from Eq. (|T2|) at the pole position: 



x 2 Tl VI 
1 - x 2 d z t n 



{a-i)i 



N- 



tN 



(15) 



leading to 



(16) 



From Eq. (fT5|) one obtains 



'TV 



S12S21 



m^r — 771,2 



2 vi 
X En 



-22 



(20) 



with the E evaluated at z = mjy- Note that for E12 = 
E 2 2 = one recovers Eq. (TTT1) . 

Again, as in the one-resonance case, we calculate x in 
order to determine the bare nucleon coupling The 
physical residue, given once more by Eq. (|13[) . has to 
agree with the residue of T p from Eq. ([5]) , 



The bare mass m b N and coupling /jy of the nucleon are 
then calculated by inserting Eq. (fTB")) into (JHJ and (TTTT) . 
which leads for to 



'A 1 



niiv - x 2 En(z ). 



(17) 



This matching procedure is valid for the nucleon in 
case it couples only to the ttN channel. In Ref. [62[, the 
renormalization for one bare s-channel state coupling to 
more than one channel is developed. In contrast, in the 
present study a second bare s-channel state is included 
that is allowed to couple to multiple channels, while the 
bare s-channel nucleon state is allowed to couple only to 
the irN channel. The corresponding formalism is devel- 
oped below. We have also developed the formalism for 
two bare s-channel states with both of them coupling to 
multiple channels. The resulting set of equations needs 
to be solved numerically for m b N and fy, but we see no 
need here to couple the nucleon pole at z = 939 MeV 
to channels with a significantly higher threshold energy 
than ttN; at the nucleon pole, the effect of omitting the 
coupling to higher-lying channels is fully absorbed in the 
value of x. 

The renormalization for the case of two bare s-channel 
states is derived in the following. The nucleon is allowed 
to couple to the nN channel only, while the second state 
couples to all channels. That second state leads to a 
resonance at z ~ 1650 MeV in the P\\ partial wave as 
discussed in Sec. IIVDI 

The nucleon pole position is given by a zero of the 
determinant of D in Eq. (JSJ) , 

(z - m\ - En)(z - m\ - E 22 ) - E 12 E 2 i = , (18) 



with zq = mjv and E, 



In Eq. ([IS]), 



m\ = m N and m 2 is the bare mass of the second genuine 
resonance; m\ is a free parameter of our model which is 
determined in the fitting process. In addition to Eqs. (flU)) 
and pip , we define the following self-energies: 



E12 = /, 

2j21 



JV ^12 d — x fnNN E 1 ^ — X Ei 2 



red 



2 ^21 



pred 
NN ^21 



X E? 



(19) 



Note that E12 (= E21) contains all bare couplings of the 
second bare state, that are, like m^, free fit parameters. 



(a-i)n 



N^ttN 



(r?, r§) 



-22 



E21 



E12 



Ei; 



771 V 



Eii)(* 



E22) — E12E21] 



(21) 



with both sides of the equation evaluated at z = m^. 
Inserting Eqs. (|10[) and (|19[) . we arrive at an equation for 
x 2 that only depends on known or fitted quantities: 



1 



G 2 E 12 (&E 22 -l) + 2d 2 Ei 2 



9 z En + G2E12 



G 2 - i r? + r^Ei2j [G^n + v^ 



(22) 



with G^ 1 = 7777V — '772 — E22- Note that for E12 = 0, one 
recovers the one- resonance case of Eq. (fTBl . 

The bare mass and coupling of the nucleon are calcu- 
lated by inserting Eq. (J22J in Eq. and Eq. $TQ) . The 
entire calculation is performed for each step in the fitting 
process. 

In the scheme derived here we ensure the correct prop- 
erties of the nucleon in the presence of a second genuine 
state. In principle, one could introduce further bound- 
ary conditions, such that the residues of the nucleon pole 
into the channels r/N and KY satisfy the same SU(3) 
constraints as the meson-meson-baryon vertices in the t- 
and it-channel exchange diagrams (cf. Appendix IB 21) . 
Those effects are, however, very small, as these channels 
are separated from the nucleon pole by several hundreds 
of MeV. Similarly, for internal consistency one could re- 
place the effective meson-meson-baryon vertices in the 
t- and u-channel diagrams by bare ones and then renor- 
malize them consistently with the nucleon renormaliza- 
tion procedure. As the t- and 7i-channcl diagrams enter 
T NP which is subsequently used for the renormalization 
of the nucleon, such an extension of the scheme would 
obviously require a self-consistent treatment beyond the 
scope of the present work. 

At this point, we would like to emphasize that bare nu- 
cleon and resonance masses and coupling constants have 
no physical meaning, because they depend on the scheme 
(in our case, the cut-off values of the form factors and 
their parameterization). 
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In addition, some bare resonance parameters tend to 
be strongly correlated: e.g., for the high- lying pN chan- 
nel, the pN self energy £ [cf. Eq. ((4])] has a small imag- 
inary part at the considered energies and thus the bare 
pN coupling to a resonance is almost fully correlated with 
the bare (real) mass. We quantitatively contrast bare and 
physical parameters, the latter being pole positions and 
residues, in Sec. HVEl and llVFl 

For the sake of reproducibility of the results, the bare 
resonance parameters are quoted in Tables IVIII and IVIIII 
in Appendix|21 Because of the reasons stated above, they 
should not be compared to lattice or Dyson-Schwinger 
results. 



III. RESULTS 
A. Data base 

The data that enter our calculation are displayed in 
Figs. I51IM1 together with our results of fits A and B of the 
analysis as explained in Sec. MI B[ Data sets which differ 
only by a few MeV in scattering energy are depicted in 
the same plot, whereas in the fitting procedure, of course, 
the exact energy value for every data set is taken. 

While for the reaction irN — > irN the partial waves 
from the GWU/SAID analysis [54[ are used, for the in- 
clastic channels it N —> r/N and ttN — > KY we fit directly 
to total and differential cross sections and to polarization 
observables. The bulk of the existing data for the inelas- 
tic channels was obtained in the 1960's and 70's. Though 
many experiments have been carried out at different fa- 
cilities, unfortunately, there are still energy regions where 
the data situation is not ideal, cf. for example the dis- 
cussion in Sec. MI Di on the reaction ir~p — > r/n. For the 
rjN, KA and K°T,° channels, there is a qualitative dif- 
ference between data on differential cross sections da/dil 
and polarizations P: in general, the polarization mea- 
surements exhibit much larger error bars and sometimes 
take on unphysical values greater than 1. For certain 
energies, there exist also conflicting measurements. The 
spin-rotation parameter /3 has been measured only for 
the K°A and the K + T, + final states, and only for a few 
energies. In the following sections, the data situation in 
each reaction is critically reviewed. The weights we apply 
to specific data sets in the fit are discussed in detail. 

Note that a measurement of da/dfl, P, and f3 rep- 
resents one possible full set of observables to disentan- 
gle the amplitudes model-independently for scattering 
of pseudo-scalar mesons off spin-1/2 baryons. In prac- 
tice, one needs to make model assumptions, as done 
here, to obtain sensible results, because the experimen- 
tal precision is far from being sufficient for a model- 
independent analysis. To disentangle the ttttN dynamic 
model-independently, one would need many more observ- 
ables anyway. 

The combined study of the reactions n + p — > A' + £ + , 
ir~p A' £ , and in particular ir~p A + £~ is es- 



sential to disentangle the isospin content of the ATE final 
state. For the A+S+ and A°S° final states, the pro- 
nounced forward peak shows the onset of ^-channel dom- 
inance which at energies > 3 GeV is most ec onomically 
parameterized in terms of Regge exchanges [l72l Il73| . 
In this respect, the A + S~ final state, although omit- 
ted in many analyses, is of particular interest, because 
it is sensitive to w-channcl exchanges, as i-channel pro- 
cesses would require the exchange of a charge 2 particle at 
the potential level. Due to its u-channel dominance and 
the small total cross section, the n~p A' + £~ process 
imposes strong constraints, despite the somewhat poor 
data situation. In other words, we found parameter sets 
that fit 7T+p -> A'+S + and ir~p — > A°S° somewhat bet- 
ter than in the present result, but once the data of the 
A' + E~ final state were included into the minimization 
procedure, the partial-wave content changed. See also 
Sec. MI III and Fig. 1331 where we test the consequences if 
the data of the reaction ir~p —> A' + £~ are ignored in the 
fit. 



B. Numerical details and fit parameters 

In the fitting procedure, we consider two different sce- 
narios. In case of fit A, we start from the parameter set 
of a preceding Julich n N model [63[ . The addition of new 
channels in the present model and the requirement that 
many more data have to be described ensure that the pa- 
rameter space is sufficiently explored in the course of the 
fitting procedure before an acceptable result is reached. 
Still, in their overall properties, the new fit A and the 
old interaction potential remain fairly similar. In par- 
ticular, the sizable renormalization of the (bare) nucleon 
mass in the old model [63[ occurs in the present fit A 
too; cf. Sec. IIVDI For fit B, we start out from a radi- 
cally different scenario so that the parameter space has 
to be searched again completely anew. The prime rea- 
son for such a strategy is that we would like to find out 
how far our results and, specifically, the extracted reso- 
nance parameters are sensitive to the starting conditions 
of the fitting procedure. Certainly, for a solid assessment 
of the uncertainty of the extracted amplitudes and pole 
properties, one would need to carry out more fits under 
different conditions. The high numerical effort, however, 
puts limits on the number of fits that can be performed. 
To make sure that the starting point is indeed quite dif- 
ferent in the fitting procedure for fit B, we imposed from 
the beginning that here the renormalization of the nu- 
cleon mass should be as small as possible. Note that the 
amount of renormalization is not an observable and thus, 
in principle, should not influence the results. It turned 
out that also for this scenario a comparable description of 
the data can be achieved. Interestingly, in this scenario 
the coupling of the ttN system to the effective ttttN chan- 
nels, notably the aN channel, had to be increased signifi- 
cantly, and the strength of the interaction within the aN 
channel increased, too. The differences in the scenarios 
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for fit A and B mentioned above imply consequences for 
the resonance generation of the Roper resonance which 
will be discussed in detail in Sec. IIVDI 

The present results were obtained in a fit procedure 
using MINUIT on the JUROPA supercomputer at the 
Forschungszcntrum Julich. To be able to handle large 
amounts of data in combination with time-consuming fit- 
ting techniques, an efficient parallclization of the code us- 
ing MPI (Message Passing Interface) was mandatory. To 
find a good local minimum, usually 50-100 parallelized 
energies were used, while for the final production runs 
300 and more processes ran in parallel. As discussed be- 
fore, the decomposition of the amplitude into a pole- and 
a non-pole part according to Eq. ^ provides the pos- 
sibility to vary resonance parameters without the need 
to recalculate T NP every time. The architecture of the 
parallelization takes advantage of this possibility to save 
CPU time: for every step in the determination of pa- 
rameters tied to T NP , all T p parameters are fully opti- 
mized, i.e. we perform a nested fit. The computation 
time can be reduced by two orders of magnitude through 
this method. 

Apart from the parallelization in energy, some pro- 
cesses perform special tasks such as the renormalization 
of the nucleon and the control of the nested fit; combin- 
ing iterative, nested fitting with the parallel architecture 
is not trivial. 

Compared to Ref. [6(|, a large amount of additional 
data for total and differential cross sections and polar- 
ization observables for the rjN and KY channels was 
considered as input to the fit. The current set amounts 
to about 6000 data points. The program is organized in 
a way that much larger amounts of data could be han- 
dled without slowing the calculation down severely. This 
will be especially important once photo- and electropro- 
duction data, that comprise more than 10 5 points, are 
included in the analysis. Also, in principle an additional 
parallelization in partial waves is possible but was not 
required so far. 

In the previous analysis [66j . the reaction ir + p — > 
K + Y, + and irN scattering were considered and only res- 
onance parameters, i.e. bare masses and couplings of the 
resonances to the different channels, were fitted. In the 
present study, in addition the important T NP parameters 
are varied. Those are the cut-offs of the form factors in t- 
and it-channel exchange diagrams, shown in Figs. [5] and 
12] and quoted in Appendix iBl 

Bare resonances with a total spin through J = 9/2 are 
included, with corresponding new parameters. One bare 
s-channel state is included in each of the isospin 1=1/2 
partial waves D13, -D15, F15, P13, Fn, H\g and Gig, while 
we have two in 5*11 and P\\. In the / = 3/2 sector, one 
bare s-channel state is included in the S31, -D33, 
P3I1 F37, G37 and G39 partial waves and two are 

included in P33. These states couple to all channels irN, 
pN , r)N , 7rA, KA and KY, if allowed by isospin. Bare 
s-channel states are included as demanded by data. The 
significance of the resulting resonances is discussed for 



each partial wave in Sec. IIVB1 

In total, we have 196 free parameters, of which 128 
are resonance parameters and 68 belong to the T NP part 
(64 cut-offs in the form factors of the t- and u-channcl 
exchanges plus 4 coupling constants). Note that varying 
a parameter that enters V NF changes the amplitude in 
all partial waves — and also for all reactions. While most 
of the T NP parameters are responsible for the part of 
the amplitude that slowly varies with energy, the Roper 
resonance is dynamically generated in the present ap- 
proach. This means that here the T NP parameters have 
to be chosen to reproduce the resonance and simultane- 
ously provide the "background" in all other partial waves. 
The values of the resonance parameters are quoted in Ta- 
bles IVIII and IVIIII and the V NP parameters can be found 
in Tables IxITlI and IXTVl 

In the course of fitting, we assigned special weights to 
some data sets to ensure that the fit respects certain fea- 
tures in the data. For example, there is a pronounced 
but narrow forward peak in the reaction ir~p — > A" A 
(cf. Fig. \W\i . At higher energies, z > 2 GcV, the peak 
is consistently present in different experiments, but al- 
ways with low statistical weight. Without giving those 
data special weight in the minimization procedure, some 
fit results would not reproduce the forward peak. Also, 
varying the data weights in the course of fitting helps 
prevent the fit from getting stuck in shallow minima of 
the parameter space. 

We do not claim a sharp upper limit in energy for the 
validity of the current approach; for example, the higher 
partial waves in ttN scattering have clear resonances be- 
yond 2 GeV [54| and in such cases it is important to 
fit even beyond the resonance position. In any case, in 
the present Lagrangian-based framework, the amplitude 
allows, in general, for a well-behaved extrapolation to 
higher energies; by contrast, in analyses in which the 
potential is parameterized purely phenomenologically in 
terms of polynomials, there may be little control on the 
amplitude outside the fitted energy region. As expected 
from the physics-driven parameterization of our ampli- 
tude, the overall agreement of the fits with the data be- 
yond 2 GeV is, for most observables, quite good. 

In the following, we discuss the data situation and our 
description of the data for the different reactions. For the 
reaction 7rA — > nN, we show the partial waves of our fit 
result. For the inelastic reactions ttN — > r]N, KY, there 
is an overall, energy- dependent phase that cannot be de- 
termined directly from experiments, and therefore we re- 
frain from giving the corresponding amplitudes. Instead, 
we show the partial cross sections, i.e. cross sections for 
each partial wave (amplitudes can be provided upon re- 
quest). Due to this unknown energy-dependent phase, 
other results should be compared to the present ones at 
the level of partial cross sections. We stress that for any 
meaningful comparison on the level of amplitudes, it is 
mandatory to take the mentioned phase ambiguity into 
account. 
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C. ttN nN 

The energy-dependent partial-wave solution of the 
GWU/SAID analysis [54[ up to H waves is used as input 
for the elastic ttN channel, both for isospin 1/2 and 3/2. 
The points displayed in Figs. [S][5] and in Appendix [Cl 
however, show the energy-independent (single-energy) 
solution of Ref. [54j . The corresponding error bars repre- 
sent uncertainties of the theoretical analysis and cannot 
be used for calculating the \ 2 - Therefore, to this point 
no true statistical uncertainties on extracted quantities 
(poles and residues) can be determined in our analysis. 
The inclusion of the actual ir N ttN experimental data 
in the future will make this possible. Note that only 
the GWU/SAID group fits to ttN -> nN data while 
in other analyses (present app roach, EB AC IjH W^L Fz2l — 
[Z|, Bonn-Gatchina [HI HH, Giefien [l36Ml39| rKent 
State (l4l fl49j . ...) their result is used as input. For 
the 20 partial waves of irN — > irN considered, we fit to 
synthetic data generated from the energy-dependent so- 
lution [54j in 5 MeV steps in scattering energy z and as- 
sign an error of 0.01 to the real and to the imaginary part. 
Although we refer to the 2006 paper [54[ throughout this 
study, we fit to the unpublished WI08 solution of the 
SAID group accessible through a web interface 54L In 



addition, we fit to the ir N scatt ering lengths of Ref. jl74j | 
and scattering volumes of Ref. |l75| . 

Figures [U to [5] show the / = 1/2 and I = 3/2 elas- 
tic ttN — > ttN partial-wave amplitudes. The alternative 
representation in terms of phase shifts and inelasticities 
can be found in Appendix [Cl The result of fit A of this 
study is represented by the red solid lines, while the blue 
dashed lines denote the result of fit B. 

The quality of the result has to be seen in perspec- 
tive: As Fig. [2] shows, there are only four tunable pa- 
rameters for the ttN — > ttN non-pole transitions, namely 
the form factor of the irNN vertex with TV in the it- 
channel, one form factor for the a- and one for the p- 
exchange, and one for the 7riVA vertex with A in the 
it-channel. Of these four parameters, two appear also in 
different transitions, e.g. irN — > ttA with A or A in the 
it-channel. Thus, these two parameters are also responsi- 
ble for the inelasticities. An additional constraint comes 
from the very small isoscalar scattering length. There is, 
therefore, only very little freedom in the fit of the non- 
resonant ttN — > ttN transition potentials. Some limited 
additional freedom comes from coupled-channel effects, 
through which also other parameters have influence on 
the non-pole ir N — > irN amplitudes. In any case, the fact 
that the non-resonant part of the amplitude in 20 par- 
tial waves can be described with very few (2-4) directly 
relevant parameters, should be considered a remarkable 
success. 

In some partial waves, there are shortcomings. E.g., for 
(Fig. \7} we achieve a good description of the imagi- 
nary part in both fits A and B. Still, the description of 
the real part is not satisfactory at higher energies. As 
Fig. [17] shows, the phase shift in S31 is well described, 
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FIG. 5: Reaction ttN -> tvN, isospin / = 1/2, S- to D- 
waves. Solid (red) lines: fit A; dashed (blue) lines: fit B; 
points: GWU/SAID partial-wave analysis (single-energy so- 
lution) from Ref. [5ll. 



and shortcomings originate from the inelasticities. Sim- 
ilarly, the description of the Su partial wave at higher 
energies could be improved. This might be due to chan- 
nels not yet included in the approach as, e.g., uN. In 
this context, one should also mention pion- induced 4> pro- 
duction on the nucleon, to be considered in the future, 
which is closel y tied to the KY channels as shown in 
Refs. fl76l fl77| . 

For many of the 7rA partial waves, fit A and B provide 
similar descriptions with only small differences, predomi- 
nantly at higher energies. The most apparent deviations 
appear in the Fn partial wave in Fig.jS] Here, because of 
the small amplitude, the resonance position is not easy 
to be fixed. Also in the Fi§ partial wave, discrepancies 




FIG. 6: Reaction irN tyN, isospin / = 1/2, F- to H- 
waves. Solid (red) lines: fit A; dashed (blue) lines: fit B; 
points: GWU/SAID partial-wave analysis (single-energy so- 
lution) from Ref. 



FIG. 7: Reaction ttN -t tyN, isospin / = 3/2, S- to D- 
waves. Solid (red) lines: fit A; dashed (blue) lines: fit B; 
points: GWU/SAID partial-wave analysis (single-energy so- 
lution) from Ref. 



between fit A and B appear at higher energies. The Pn 
partial wave is of special interest and will be further dis- 
cussed in Sec. IIVDI 



D. Reaction nN 



r]N 



Most data of the reaction ir~p — > r\n date back to 
the 1960's to 1980's and exhibit inco nsistencies among 
each other as scrutinized in Ref. jl78l | where a detailed 
overview and analysis of the different experiments and 
the quality of the corresponding data is given. See also 
the selection of rjN data provided by the GWU/SAID 
group [53l | and in Ref. [l79j . Due to a frequent un- 



derestimation of systematic errors and issues concerning 
beam calibration and normalization, a careful selection 
and rating of the data sample is necessar y. W e mainly 
follow the esti mations made thereto in Ref. [178| , see also 
Refs. [73l.fl50|. 

In the pr esent an alysis , differential cross section data 
from Refs. [180| | to 181 for 38 different energies from 
z = 1489 MeV up to z = 2235 MeV are included in 
the fit . Fo r the more recent exp erim ents by Prakhov 
et al. |l 80j ] and Bayadilov et al. [l82j . we employ the 
syste mat ic uncertainties given by the authors. Following 
Ref. |l83| . a systematic error of S sys = 10% of the exper- 
imental value is add ed to the statistical error of the data 
from Morrison et al. [l84| . As Fig.[9]shows, the Morrison 
data (empty circles) are in conflict with the Prakhov data 
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FIG. 8: Reaction nN -» tyN, isospin I = 3/2, F- to if- 
waves. Solid (red) lines: fit A; dashed (blue) lines: fit B; 
points: GWU/SAID partial-wave analysis (single-energy so- 
lution) from Ref.[5l|. 



(filled squares). We have decided to rely on the Prakhov 
data in this analysis and drastically reduced the weight 
of other data in the fitting procedure, in case they are 
in conflict. E.g., for the data sets at z = 1498 MeV and 
z = 1499 MeV, the data by Prakhov were weighted with 
a factor 25 compared to Morrison in both fits. 

The uncertainties stated in Kozlenko et al. [l85j are 
predominantly of statistical nature. Therefore, a system- 
atic error of S sy s = 10% was assigned. For Debenha m et 
al. [l86j we use the same uncertainties as Refs. (zj| Il50| 
and employ a systematic error of 10% and 0.02 mb. A 
systematic error of S sys = 11.4% of the experimen- 
tal value was add ed to the statistical error of the data 
from Deinet et al. [l87j | as advised by the authors them- 
selves. The large energy interval from z = 1636 MeV 



to z = 2235 MeV is almost solely covered by Bro wn et 
al. 181 1 , supplemented by data from Richards et al. |l88j 
at medium energies. For Richards et al., we assume an 
additional error of 8 sys = 9% to 14% as given b y th e 
authors in the original paper. As argued in Ref. [178|, 



the data from Brown et al. [I8lj are highly question- 
able at low energies due to a miscalibration of the beam 
momentum and, therefore, are not used as input in the 
present calculation. At higher energies, on the other 
hand, the data might be used if the beam momenta are 
lowered by 4%, together with an additional systematic 
error of 10%. In any case, the Brown data enter the 
fit with a much lower weight and we make no attempt 
to fit tho se data quantitatively (in the recent analysis of 
Ref. jl48l | the data are quantitatively fitted) . For the to- 
tal cross section, we include the data that are rated as 
reliable in the GWU/SAID analysis [H, H3|. These are 
shown with filled squares in Fig. [TTJ Note that in the 
GWU/S AID approach, the total cross-section data from 
Prakhov [180| are not included to avoid double counting 
with the respe ctive differential cross-section data that are 
included jl89l |. 

The results for the differential cross section are shown 
in Figs. M and [TU1 where the red solid line denotes fit 
A and the blue dashed line represents fit B. The over- 
all agreement is good for both fits and only slight dif- 
ferences in the data description can be observed. For 
the lowest energies, it should be noted that due to the 
non-mon ochr omatic pion beam and uncertainties in cal- 
ibration |l80| , in combination with the very steep rise of 
the cross section, the errors on the data are la rger than 
those shown in Fig. [HI see the discussion in Ref. |45l . ll80| . 

For scattering energies of z ~ 1575 MeV and higher, 
the data situation becomes problematic. In Fig. [TT] the 
rjN total cross-section data reflect these problems. The 
data quality for the differential cross sections is quite low 
and conflicting at these energies, and it is difficult to re- 
liably fix the amplitude. Thus, fit A (red solid line) and 
fit B (blue dashed line) give different descriptions of the 
total cross section from z ~ 1700 MeV to z ~ 1850 MeV, 
where fit A seems to be too low. This aspect also occurs 
in the differential cross sections at the corresponding en- 
ergies, which are better matched by fit B, although the 
weight on the different data sets was comparable to the 
one in fit A. 

Polarization data for the reac tion ir~p — > rjn have only 
been published by Baker et al. jl9l| in form of 12 polar- 
ized cross sections from z = 1740 MeV to z = 2235 MeV. 
Here, in principle, the same beam calibration problems 
as in Brown et al. [18lj could occur, since the same ex- 



perimental apparatus and a similar beam line were used. 
Moreover, the results from the latter experiment enter 
the analysis by Baker in the event selection of the final 
data sample. In Fig. [T^] we still show results for the po- 
larized cross section. The form of the polarization does 
not depend strongly on the energy, so that a 4% error 
in the beam momentum would not be outright notice- 
able. Because of the problems mentioned, the polarized 
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FIG. 9: Differential cross section [1/2] of the reaction n~p — > 
r\n. Solid (red) lines: fit A; d ashed (blue) lines: fit B; data : 
filled squares from Ref. Il80|j : filled circles from Ref. [l82t | : 
empt y cir cles from Ref. Il84| : empty triangles up from 
Ref. Il85|; stars from Ref. [1861 ]: fill ed t riangles up from 
Ref. [187| ; em pty squares from Ref. [188| | ; fille d diamonds 
from Ref. 190]; empty diamonds from Ref. [I8l| . 



cross sections entered our analysis with a very low weight 
although the fits shown in Figs. [T21 and [TBI exhibit an ac- 
ceptable description of the data. 

To summarize, given the problematic data situation for 
the reaction 7r~p — > rjn above z > 1.56 GcV, we found it 
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FIG. 10: Differential cross section [2/2] of the reaction n~p — > 
r\n. Solid (red) lines: fit A; dash ed (blue) lines: fit B; 
data: em pty squares from Ref. [1881 ]: empty diamo nds from 
Ref. [l8l| . Note that the energy for data from Ref. |l8ll ] has 
been lowered (see text). The latter data are unreliable and 
enter the fit with very low weight. 



difficult to pin down the partial-wave content at higher 
energies. 

For reasons stated at the end of Sec. IIII B[ we refrain 
from showing amplitudes for the inelastic reactions. In- 
stead, we present the partial cross sections in Fig. [TH 
where the solid lines represent fit A and the dashed lines 
denote fit B. 

As can be seen in Fig. [TH the strong increase of the 
total cross section just above the r\~N threshold [z = 
1486 MeV) is dominated in both fits by the Sn wave, 
which is in line with the known strong branching ratio 
of the iV(1535) l/2~ resonance to the r/N channel (cf. 
Table H]T[) . The influence of the Sn partial wave remains 
large over the whole energy range. At z ~ 1.65 GeV, the 
Pn partial wave shows noticeable impact but decreases 
quickly. This resonance plays also a ro le in 77 photopro- 
duction, as recently shown in Ref. jl40| . In addition, the 
D15 partial wave is strong at z ~ 1.7 GeV. At higher en- 
ergies, the P13 and D13 gain some influence besides the 
Sn. Higher partial waves like the Fn, G17, Gig and Hig 
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FIG. 11: Total cross section of the reaction n~p — > r\n. Solid 
(red) lines: fit A; dashed (blue) lines: fit B; data: the filled 
squares indicate experiments accepted by the GWU/SAID 
group [S^l; open circles from Prakhov et al. [180(] . 
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FIG. 12: Polarization [1/2] of the reaction n~p — > r\n. 
Solid (red ) lin es: fit A; dashed (blue) lines: fit B; data 
from Ref . [I9ll ] . Note that the beam energies for data from 
Ref. |l9ll ] have been lowered and the data are included in the 
fit with low weight (see text). 



are very small overall. The partial-wave content of the 
total cross section is similar in fit A and B. Note the loga- 
rithmic scale of Fig.QJJ which causes discrepancies of the 
two fits for the small partial waves to appear enlarged. 
The differences in the Fyj contribution, in contrast, are 
already apparent in the elastic irN channel, c.f. Fig. [B] 



FIG. 13: Polarization [2/2] of the reaction ir~p — > r\n. Solid 
(red) line s: fit A; dashed (blue) lines: fit B; d ata from 
Ref. [l9ll ]. Note that the energy for data from Ref. |l9l| has 
been lowered, cf. Sec lIII Al The data for P were not included 
in the fit. 
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FIG. 14: Partial cross sections of the reaction ?r p -> r\n. 
Solid lines: fit A; dashed lines: fit B; data: as in Fig. 1111 



The u-shape of the differential cross section visible from 
threshold to around z ~ 1.6 GeV is due to the SD-w&ve 
interference, as can be seen in Fig. [T5l 



E. Reaction 7r 



P 



K°A 



For the reaction ir~p —> K a A, total and differential 
cross sections, polarization, and spin rotation parameter 
measurements from z = 1626 MeV up to z = 2405 MeV 
are included in the fit. The data base contains differential 
cross sectio ns at 46 energ ies and 8 different experiments 
from Refs. |l09l . Il92l4l98| . In contrast to the r/N case, 
no severe inconsistencies appear in the KA differential 
cross sections and systematic uncertainties were applied 
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FIG. 15: Differential cross section of the reaction f p-> r\n 
at z = 1512 and z = 1534 MeV. Only results of fit A are 
shown. Solid (red) lines: full result; dashed (green) lines: Su, 
Pn and P13 contributions only; dash-dotted (black) lines: Su 
and D13 contributions only; data: as in Fig. [9] 

as given by the respective authors in the original publi- 
cations. 

For the recoil polarizatio n, we use the data of 
Refs. PM EM [Hi [HI [HI- 27 energies are covered. 
At backward angles > 90°, the data exhibit incon- 
sistencies and huge error bars. It should be noted that 
several data points take unphysical values larger than 1. 
A measuremen t of the spin-rotation parameter ft was re- 
ported in Ref. 199]. All seven data sets at energies from 
z = 1852 MeV to 2262 MeV are included in our analysis; 
cf. Fig.[T9l Note that this observable is 27r cyclic. Errors 
were applied as stated by the authors. The spin-rotation 
angle can also be written in terms of the spin-transfer 
coefficients -Ky, 

_ 21m (gh*) _ K zx _ K xz _ K zx _ K xz 
~ M 2 - \h\ 2 ~ K zz ~ K xx ~ K xx ~ K zz ' 

(23) 

where g is the spin-non-flip and h the spin-flip amplitude 
(see Ref. [66| for a definition of observables and ampli- 
tudes). 

The data for total cross secti on of KA production in - 
cludcd in the fit are from Refs. [109L 
Additionally i nclu ded data from Ph.D theses are refer- 
enced in Ref. [205j . 

The current fit result is compared with the data in 
Figs. QH to [201 Tlie description of the n~p -> KA differ- 
ential cross section, polarization, and spin rotation pa- 
rameter is satisfactory for both fits A and B. At low en- 
ergies, the description of the differential cross sections 
are very similar in fit A and B, while at medium energies 
fit B seems to be a little better. At higher energies, the 
peak in forward direction is more pronounced in fit A, 
but also clearly visible in fit B. The backward peak, on 
the other hand, is better described in fit B. Note that the 
behavior of fit A at backward angles for higher energies, 
where it bends downwards at cos ~ — 1, also appears in 



some data sets, e.g. at z = 2159 MeV and z = 2182 MeV. 
Here, the data situation seems to be slightly inconsistent. 
Differences in fits A and B for the polarization (Fig. [T7|) 
around cos© = —0.5 show that the data quality is not 
good enough to reliably fix the amplitude. The total cross 
section, shown in Fig. [201 is a l so wei l described, although 
there are conflicting measurements around z ~ 1.7 GeV. 
This is reflected in the discrepancies around that energy 
between fits A and B. Neither of the two fits, however, 
reaches the highest data points. It should be noted that 
the differential cross section in that energy range, shown 
in Fig. [Tni is well described; this applies especially to fit 
B. Hence, our result suggests that the total cross section 
peaks with a ~ 0.7 mb, and that the conflicting data 
with higher values should be discarded. 

In Fig. [19j we show the description of the spin- rotation 
parameter ft. The most apparent discrepancies in fits 
A and B appear at backward angles where no data are 
available to constrain the fit. This demonstrates the need 
for precise data in all observables and angles to reliably 
determine the amplitude. 

The partial cross sections are shown in Fig. [2U The 
solid lines indicate fit A and the dashed lines represent 
fit B. As will be discussed in Sec. IIV Dl the description of 
the K°A data at low energies requires a Pn resonance, 
which, in turn, induces a structure in the 7rA^ — > ttN am- 
plitude that matches the single-energy solution [54j even 
though the (unstructured) energy dependent solution in 
ttN — > ttN was fitted. The strong influence of the Pn on 
the K°A channel at low energies is reflected in the par- 
tial cross sections. Of the remaining partial waves only 
the Pi 3 shows some influence at low energies, which is, 
however, not comparable to the Pn and Su. At higher 
energies above 2 GeV, the Pn, Su and Pi 3 decrease and 
the D13 and Fn become stronger in both fits. 



F. Reaction 7r~p K°YP 

The results for the reaction 7r~p — > K°T,° are shown 
in Figs. (55] to US] Here, we take into account differential 
cross sections at 29 different en ergies from z = 1694 MeV 
to z = 2405 MeV from Refs. EM E13, Hoi Hq3- 
References [206l l207l | also provide polarization data. The 
two data sets, shown in Fig. [221 do not seem to be con- 
sistent at all places and both exhibit large error bars, 
which reduces the impact on the fit result. Also, in 
several places the data take unphysical values larger 
than 1. For all data in this channel systematic uncer- 
tainties were applied as given in the original publica- 
tions. The data for the total cross s ection are taken from 
Refs. fl93l. EM E M [203 1 f205l - [209 | [ ^ We do not use the 
data provided by Livanos et al. for this reaction 

and other KH final states as they carry large error bars 
compared to other data. 

The current fit results match the differential cross sec- 
tion well. There are only minor differences in fits A and 
B, predominantly at higher energies where the forward 
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FIG. 16: Differen tial c ross section of the reaction 7r~j 

circles fro m R ef. |l95l |: empty squares from Ref . |l96l |; partially filled circles from Ref . [I92t [: stars from Ref. [19 
from Ref. |l93f | : empty diamonds from Ref. [I09f j : partially filled squares from Ref. 194]; filled diamonds from Ref. 



K°A. Solid (red) lines: fit A; dash ed (blue) lines: fit B; d ata: filled 

; tr iangles 

in. 
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FIG. 17: Polarization of the reaction ir~p — > K°A. Solid (red) 
lines: fit A ; dashed (blue) lines: fit B; dat a: fil led circles from 
Ref. [192| ]; emp ty triangles up from Ref. Il93|] ; filled squares 
from Ref. [199| ; empty diamonds from Ref. |l09| . 



peak is more pronounced in fit B. Around z ~ 1.8 GeV, 
the total cross section seems to be underestimated a bit 
by both fits, but the inspection of the differential cross 
section in Fig. [M] shows that there are only a few data 
points at backwards angles that are not matched well. In 
addition, the distribution of the data at backward angles 
at z = 1792 MeV and z = 1797 MeV and the adjacent 
energy bins suggests inconsistencies in the data. 
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FIG. 18: Polarized cross section of the reaction 1 p-> K°A. 
Solid (red ) lines: fit A; dashed (blue) lines: fit B; data from 
Ref. fl95lt. 
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FIG. 19: Spin-rotation parameter /3 of the reaction iv~p — > 
K°A. Soli d (re d) lines: fit A; dashed (blue) lines: fit B; data 
from Ref. [199^ . Note that ft is 2ir cyclic. 



Next, we discuss the partial-wave content of the reac- 
tion n~p — > K°YP. The partial cross sections a L ' 2I: 2J for 
the other reactions ir~p —> K + Y,~ and ir + p — > A' + S + , 
where L( J) is the orbital (total) angular momentum and 
I the isospin, are related to the partial cross sections of 

0v0 



the 7r p-> K"TP reaction, 



-"1;2J 

T K+J2- 

T K+Y. J ' 



9 



'A'+S- 



^3;2J 

2 



(24) 



We therefore only show the partial cross sections for the 
reaction n^p -> K°T,° in Figs. [25] and [26] The partial 
cross sections for the final states K + Y;~ and K + Y, + can 
be obtained from these curves and Eq. (j2"4"|) . 

In contrast to K°A, the isospin 7=1/2 partial cross 
sections for the KT, channels (Fig. l25j) are dominated 
almost entirely by the S\\ partial wave in both fits. The 
P\i partial wave plays nearly no role at all in fit B, while 
in fit A it is second to the influence of the S\i at very 
low energies. At medium energies, the D15 and Pi 3 are 
second to the influence of the Sn in both fits, but not of 
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FIG. 20: Total cross section of the reaction n~p — > K°A. 
Solid (red) lines: fit A; da shed (blue) lines: fit B; data : 
filled circles from Ref. |l92t |; fille d squares from Ref. [l09l |: 
empt y di amonds from Ref. |193| ]: empt y tr iangles up from 
Ref. jl94l ] ; filled trian gles up from Ref. [l95| | ; filled tri angle s 
down from Ref. [200l |; empty tria ngles down from Ref. [20ll | : 
filled tria ngles right from Ref. [l96l ]: empty t rian gles left 
from Ref. [202l |: e mpty triangles right from Ref. |203||; empty 
squares fr om R ef. [I97t ] : filled diamonds from Ref. [19 
from Ref. 204]; for empty circles and crosses see Ref. 2051 ] 
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FIG. 21: Partial cross section of the reaction ty p —} K°A. 
Solid lines: fit A; dashed lines: fit B; data: same as in Fig. 1201 



comparable degree. At higher energies, they exceed the 
5x1. Note that Fig. 1251 suggests that the -D15 and Pi 3 
partial waves have the same threshold behavior for fit 
A, which is, however, not the case as a closer inspection 
shows. 

For the isospin 1 = 3/2 sector (Fig. the P33 is the 
dominating partial wave over almost the entire energy 
range up to z = 2.2 GeV, except for a narrow region just 
above threshold where the S31 is stronger than the P33. 
This distribution applies to both fits. In the present anal- 



FIG. 22: Total cross section of the reaction ir~p — > K°T,°. 
Solid (red) line: fit A; d ashed (blue) line: fit B; data: e mpt y 
triangles up: Ref. [193H : empty triangles right: Ref. |194[ |: 
filled squares: Ref. l203t ] : empty triangles l eft: Ref. [l98| ; filled 
trian gles up: Ref. [2081 ] filled circles : Re f. [2061 ]: empty circles: 
Ref. [2091 ]: filled triangles left: Ref. [207t ]. For empt y tria ngles 
down, filled triangles right, and empty squares see [2051 ]. 
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FIG. 23: Polarization of the reaction tt'p -» A'°E°. Solid 
(red) line s: fit A; dashed (blue) lines : fit B; data: circles from 
Ref. HHl; diamonds from Ref. [2071. 
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FIG. 24: Differ ential cross section of th e rea ction ?r p-> K°T,°. Soli d (red) lines: fit A; da shed (blue) lines: fit B; da ta: circles 
from Ref. |206l]: up triangles from Ref. [l93l |: squares from Ref. [l94f | : diamonds from Ref. [2071 ] : stars from Ref. [l97l |. 
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FIG. 25: Partial cross sections of the reaction ir^p — > K°Y,° 
in the isospin 1/2 sector. Solid lines: fit A; dashed lines: 
fit B. Note that the partial cross sections for the other KT, 
final states can be obtained from these curves, Fig. 1261 and 
Eq. dHJl. 



FIG. 26: Partial cross sections of the reaction %~p — > K E 
in the isospin 3/2 sector. Solid lines: fit A; dashed lines: 
fit B. Note that the partial cross sections for the other KY, 
final states can be obtained from these curves, Fig. 1251 and 
Eq. ((241). 
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ysis, around z = 1950 MeV the F37 partial cross section 
increases to the same magnitude and becomes even larger 
than the 531 . In the previous analysis of Ref. [66| , we also 
found an increase of the F37 at the same energy. The in- 
fluence of the remaining partial waves on the total cross 
section seems to be limited, with the exception of the D35 
which shows a similar development as the F37, but with 
a weaker slope. The isospin 7 = 3/2 partial- wave content 
is very similar for fits A and B. The differences become 
larger only at very high energies. Note the logarithmic 
scale of Figs. [25] and [26j which causes the discrepancies 
between the small partial waves to appear enlarged. 



G. Reaction n p 



As mentioned in Sec. IIII Al this reaction is of particu- 
lar interest due to its u-channel sensitivity. Considering 
the A + S~ final state is also crucial to reliably disentan- 
gle the isospin content of the AS final state. See also 
Sec. IIII HI and Fig. 1331 where we test the consequences if 
the data of the reaction ir~p — > A + S _ are ignored in the 
fit. 

The data base for the reaction — > A + S~ com- 
pris es only 15 differential cross section sets jl94l . ll98l . l21ll 
12121 ] from z = 1739 MeV up to 2405 MeV. Although no 
severe inconsistencies are observed, the influence of this 
reaction in the fitting procedure is limited due to the 
small number of data points. For the total cross sec tion, 
data from Refs. pH, HH, Hq! Ull Hq! ImU^ are 
employed. The data are shown in Figs. [27] and [251 The 
data base is quite limited, yet the fact that the total cross 
section is significantly smaller than for the A + S + final 
state makes the description of this channel challenging, 
setting tight limits on the size of u-channcl exchanges, 
and through rescattering, also of t-channcl contributions. 

The present description is good and very similar in fits 
A and B at low energies but shows deviations at higher 
energies, as visible in Fig. [27] For intermediate energies 
around z ~ 1.9 GeV, the backward peak from u-channel 
exchanges is underestimated, in fit B even a little more 
than in fit A, while at higher energies the present results 
tend to overestimate the data at very backward angles. 
This is also visible in the description of the total cross 
section in Fig [28j The data could be conflicting around 
z ~ 1.9—2.0 GeV and the current result clearly represents 
a compromise in the description of the data at medium 
and high energies. 

As explained in Sec. IIII FL the partial-wave content of 
the reaction ir~p — > A' + S~ can be obtained from the one 
of ir~p —> A^°S° (Figs. [231 12"6| . scaled by different isospin 
factors: as Eq. ([24"]) shows, the A + S~ final state is more 
sensitive to the isospin 1 = 1/2 sector than A°S°. This 
is important for the determination of the AS branching 
ratios of / = 1/2 resonances. 

As this reaction channel is so important, better data 
would mean a significant step forward to disentangle the 
reaction dynamics in AS production. This is further 
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FIG. 27: Differential cross section of the reaction tt p — > 
K + Y,~ . Solid (red) lines: fit A; dashed (blue) lines: fit B; 
data: circles from Re f. [21 1| : up triangles from Ref. |212t |: 
squares from Ref. [l94f | : diamonds from Ref. [l98j . 
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FIG. 28: Total cross section of the reaction n~p — > K + Y,~ 
Solid (red) line: fit A : das hed (blue) line: fit B; da ta: empty 
triangles down: Ref. |l94l] : filled circles: Ref. [203l l; filled tri- 
angles down: Ref. Il98l| : filled triangles left: Ref. 



2091 ]: empty 

triangles up: Ref . [213l |: filled squares: Ref. [21 If : filled dia- 
monds: Ref. l212ll. For filled triangles up and filled triangles 
right see Ref. [205l |. 
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motivated by the f act that high prec i sion KT, photopro - 
duction data exist flToL [TlTl IllMTM llTMl^lml . [l23l - 
Il27l ] and a global analysis of pion- and photon-induced 
reactions would greatly benefit from better data in the 
ir~p — > K + T~ reaction. Efforts in this direc tion are be- 
ing made at the E19 experiment at J-PARC [l 34j | . 



H. Reaction n+p K + T, + 

For the reaction n + p — > K + T + , the same data sample 
as in our previous analysis [66j is used, which contains 
differential cross sections and polarization measurements 
in a range from z = 1729 MeV to 2318 McV for 32 
diff erent energies. We employ the data by Candlin et 
al. 104| , supplem ented by additional data at lower en- 
ergies from Refs. [214l4218| . The latter are compatible 



with the data from Ref. jloJ ] in the overlapping energy 
region but usually have larger error bars. For higher 
energies, the polarization was re-measured in Ref. 



and found to be consistent with that of Ref. [104j up 
to small deviations. In addition, value s for the spin- 
rotation parameter (3 by Candlin et al. [106] arc avail- 
able at two energies and included in the fit. During the 
measurement of the spin-rotation parameter, the polar- 
ization was re-measured and consistency with the result 
of Ref. |104| wa s found once more . Total c ross section 

data from Refs. fio i [20 I I2H - I2H [220II225I are also in- 
cluded in the fit. No additional systematic errors beyond 
the original ones were assumed for this channel. 

The present results are shown in Figs. UHl to [M] with 
the solid (red) lines representing fit A and the dashed 
(blue) lines representing fit B. Even though we add a 
large amount of new data and extend the fit to four more 
channels, we achieve an equally good result for the dif- 
ferential cross sections, polarization and spin-rotation pa- 
rameter in the i: + p — > K + T + channel as in our previous 
analysis [6(| . Fits A and B show only very small discrep- 
ancies. A slight systematic deviation becomes visible in 
the total cross section around z ~ 1.9 GeV. The inspec- 
tion of the differential cross section around that energy, 
shown in Fig. [2H reveals that the underestimation of the 
data is mostly due to the pronounced forward peak in 
the data. In both of the current fits this forward peak is 
underestimated at this energy, but well described beyond 
2 GeV. 

As a test, we have removed the observables of the 
K + T~ reaction from the fitting procedure. This would 
correspond to the situation in many analyses where the 
data of the reaction n~p —> K + T~ are not considered. 
With this change, a refit of T p starting with fit B is 
performed, and the outcome for the total cross sections 
of the KT final states is shown in Fig. [35] (dash-dotted 
lines vs. the original fit B shown with the dashed lines). 
As the fit of the reduced data base shows, the K + T~ 
final state is, of course, not well described any more. In 
contrast, it becomes now possible to describe the maxi- 
mum of the K + T + cross section. This finding emphasizes 




FIG. 29: Differential cross section [1/2] of the reaction ir + p — > 
K+E+. Solid (red) lines: fit A.; dashed (blue) lines: fit B; 
data: Ref. [Toil, except: z = 1729, 1757, 1789 MeV from 
Ref. [Hi]; z = 1732, 1783, 1813 MeV from Ref. l215ll: z = 
1764 MeV from Ref. [HI; z = 1790 MeV from Ref. I217H. 



again the importance to include the data of the reaction 
■n-p K + T~. 

The partial-wave content of the reaction ir + p — > 
K + T + can be read off Fig. [251 on ly scaled by a dif- 
ferent isospin factor [cf. Eq. ((24)) ]. These results su- 
persede those of Ref. |(5(|. The P33 partial wave domi- 
nates starting almost from threshold up to z ~ 2.2 GeV. 
The F37 partial wave with the A(1950) 7/2 + is also well- 
determined and important while other partial waves start 
to become relevant only at higher energies. 



I. Scattering lengths and volumes 

Pion-nucleon scattering lengths and volumes of the 
current fit result arc shown in Tabic [I] and compared 
to values in the literature. The notation is , where 
L is the orbital angular momentum, and an upper plus 
(minus) sign indicates the isoscalar (isovector) scatter- 
ing length/volume. The lower ± sign indicates the cases 
J = L ± S where J (S) is the total angular momen- 
tum (nucleon spin) . Of particular interest is the isoscalar 
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FIG. 30: Differential cross section [2/2] of the reaction ir + p — > 
K + Y, + . S olid ( red) lines: fit A; dashed (blue) lines: fit B; data 
from Ref. [jo5 j. 
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FIG. 31: Total cross section of the reaction n + p — > K + ~E + . 
Solid (red) line: fit A; d ashed (blue) line: fit B; data: filled 
triangles rig ht: R ef. [l04j ] : filled squares: Re f. |214ll : filled dia- 
monds: Ref. l215fl; empty triangles up: Ref. [22y| ; filled trian- 
gles down: Ref. [2ld | ; emp ty triangles down: Ref. [22l[ |: empty 
triangles left: Re f. [222l ]: empt y cir cles: Ref. [223l |: empty 
squares: Ref. [2241 ] ; crosses: Ref. [225l | . For filled circ les, filled 
triangles left and empty triangles right see Ref. [205l |. 



TABLE I: Scattering lengths and volumes a~i± in units 
of 10~ 2 / m n an d 10 -2 /m^., compared to Ch PT calculations 
from Ref. [l74| (scattering lengths) and Ref. [l75l ] (scatt ering 
volumes). Note that for the scattering volumes in Ref. [175[ | 
no uncertainties but three different fit results are given. Here 
we quote fit 1 together with the maximal deviation among 
the three fits. 



scattering length Qq + that was determin ed to high accu- 



racy recently in the calculation of Ref. [174J ] using chi- 
ral perturbation theory. Although in the pr esent ap- 
proach we cannot match the value from Ref. [174| . we 
can indeed achieve a very small value for the isoscalar 
scattering length in both fits. This is not trivial, given 
that we have only four parameters for the non-resonant 
■nN — > TtN transitions. Those are the only parameters 
to directly tune the irN background in 20 partial waves 
from ttN threshold to beyond 2 GeV in scattering energy 
(cf. discussion in Sec. IIIIC[) . Moreover, in contrast to 
the calculation of Ref. [174| , we do not take into account 
mass differences between tv^ and 7r°, and also no elec- 
tromagnetic ef fects , and thus cannot exactly match the 
values of Ref. jl74f . For the isovector scattering length, 
our val ues a re slightly higher than the ones calculated 
in Ref. [l74{ . The scattering volumes are quite well de- 
scribed com pared to the fourth-order ChPT calculation 
of Ref. " 



175 



In Table [TT1 we show the scattering lengths of other 
channels and compare the irN and rjN values to those 
calculated in Ref. [451 ] within a unitary extension of chiral 
perturbation theory. Fits A and B from our analysis yield 
similar results. While the agreement between the results 
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11.88 
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-7.34 


-7.87 


-7.99 ± 0.335 



of fit A and of Ref. [45[ is reasonable in case of the ttN 
scattering length, for the real part of the r/N scattering 
length discrepancies can be observed between the current 
results and those of Ref. [45[ . 

The present values for the r/N scattering length can 
be also c ompa red to the compilation of results given in 
Refs. [53l. Il79j ]. The real parts of different analyses usu- 
ally differ substantially while the imaginary parts are bet- 
ter determined. The present results for the latter are in 
line with other determinations. 
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TABLE II: Scattering lengths in (in fm) of the present approach (upper rows), compared with the ChPT calculation from 
Ref. [45| (bottom row). 
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FIG. 32: Polarization [1/2] of the reaction tt+p -> K+T.+ . 
Solid (re d) lines: fit A; dashed (blue) lines: fit B; data : 
Ref. [To3], except: z = 1729, 1757 1789 MeV from Ref. [214T |: 
z = 1783, 1813 MeV from Ref. [HI]; z = 1764 MeV from 
Ref. [HI; z = 1790 MeV from Ref. [H3]; z = 1732 MeV 
from Ref. En!. 



IV. RESONANCE SPECTRUM 



A. Resonance properties 




In the literature [71|, |169{ , the decomposition according 
to E q. pi has been used since long ago. However, in 
Ref. [64| it is shown that this decomposition into pole 
and non-pole part is problematic in the sense that the 



FIG. 33: Polarization [2/2] of the reaction n + p -> K+T.+ . 
Solid (re d) lines: fit A; dashed (blue) lines: fit B; data: 
Ref. Hoi. 



pole part contains finite constant terms ~ [z — zq)° at 
the pole position z = zq. What is well defined is the 
expansion of the fully dressed, physical amplitude in a 
Laurent expansion around the pole, with the residue term 
containing all contributions of the pole, and all higher 
order terms accounting for the background. 

Thus, a resonance is uniquely characterized by its 
pole position in the complex energy (z) plane, the pole 
residues for the various channels, and the Riemann sheet 
on which it is located. 

Poles can appear on different Riemann sheets (except 
for the physical sheet of the lowest-lying channel). For 
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FIG. 34: Spin-rotation parameter (3 of the reaction n p — > 
K+T.+ . Solid (red) lines: fit A; dashed (blue) lines: fit B; 
data: Ref. [io^ |. Note that /3 is 2tt cyclic. 
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FIG. 35: Refit (black dash-dotted lines) of a reduced data 
base after taking out the observables of the reaction n~p — > 
K+T,- (right) and giving more weig ht to the K+Y.+ final 
state (left). The dashed (blue) line shows fit B, for which the 
ir~p — > K + Yi~ observables are included. 
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example, a resonance situated above the ?7-/V threshold, 
with a sizable coupling to the ttN channel and vanishing 
coupling to the r]N channel would correspond to poles 
on both the physical and the unphysical sheet of the 
77./V channel, at approximately the same pole position Zq 
(while both poles would be on the unphysical 7riV sheet). 
However, not all poles on different sheets bear physical 
interest but usually only the one on the sheet closest to 
the physical axis. For an example where a pole on a more 
distant sheet plays a role, one can consider the case of the 
iV(1535) l/2~ resonance and its interplay of poles lead- 
ing to the cusp at the r]N threshold [36|, [65[ . Such poles 
arc sometimes referred to as shadow poles. The selection 
of the part of the Riemann sheets which is close to the 
physical axis can be achieved by rotating the right-hand 
cuts of all channels towards the negative Im z direction, 
which is the choice adopted in this study. 

Apart from threshold openings on the physical axis 
at z = mi + mi for the channels with stable particles 
with masses m\ and mi , i.e. 7riV, r]N, KA, A'S, there 
are threshold openings at z = 2m 7r + mjy for the effec- 
tive 7T7rAf channels 7rA, crN, pN that also induce branch 
points in the complex plane at z = m s ± z r , where m s 
is the mass of the stable spectator particle and z r is the 
complex resonance pole position of the unstable particle. 



FIG. 36: Pole positions for the 1 = 1/2 (above) and 1 = 3/2 
(below) resonances. Squares: fit A; diamonds: fit B. Also, 
the branch points of the amplitude are shown (green crosses) . 
Note that all cuts, starting at the branch points, are chosen 
in the negative Im z direction, which uniquely defines the 
Riemann sheet of the displayed poles. 



For the cuts associated with the 7r7rA^ channels, we choose 
the same convention as for the channels with stable par- 
ticles, i.e. the cuts are rotated towards the negative Im z 
direction. An overview of the various branch points can 
be found in Fig. [551 an d Fig. [35] shows the position of 
the cuts. More details of the analytic structure can be 
found in Refs. [6^, l89j |. The resonance properties have 
been determined using the methods of analytic continu- 
ation developed in Ref. [65| and residues are calculated 
following the procedure in Appendix C of Ref. (66j . 

The choice of cuts defines uniquely that part of the 
Riemann sheets where all poles extracted in this study 
are situated (unless stated otherwise). We refer to this 
part of the Riemann sheets as second sheet. The strength 
of a pole with respect to the different channels is charac- 
terized by the residue a_i jA1!> as it appears in the Laurent 
expansion of the transition amplitude T^J from channel 
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(i to v on the second sheet, 



+ + 0{z - z ) 



(25) 



For n channels, since residues factorize in channel space, 
one has n 2 different residues, 



(26) 



with a unique set of n parameters g^ up to one undeter- 
mined global sign . To connect to the quantities quoted 
in the PDG [2261 ] . the dimcnsionlcss meson-baryon am- 
plitude r is expressed in terms of the T-matrix in the 
current normalization, 



in/ ■ 



(27) 



where k fi (E^,Lu tl ) is the on-shcll three- momentum 
(baryon energy, meson energy) of the initial or fi- 
nal meson-ba ryon system p,. The so-called normalized 
residue (NR) |226j is defined for t^, i.e. with Eq. (|2"7) 
we obtain for the transition irN — > p 



{NR)„ 



N- 



(r to t/2) 



(28) 



with the resonance width r to t = — 2Imzo [the additional 
minus sign with respect to Eq. (|27[) comes from the con- 
vention of Eq. (|3"Tj) below] . 

Defining the branc hing ratio into channel p in the same 
way as by the PDG [226j . 



l^iV-wrJV | (rtot/2) ' 



(29) 



where r is the residue of r(z), it is straightforward to 
show that the transition branching ratio equals the mod- 
ulus of the normalized residue, 



r l/2 p l/2 

= \(NR)„ N . 

t tot 



(30) 



In this study, we quote the normalized residues for ttN — > 
p instead of branching ratios, because individual branch- 
ing ratios into channels p (other than irN) are less re- 
liably determined than the quantities for the transition 
reactions irN — > p. 

As for the countin g of the residue phase 8 j226| | , we use 
Holder's convention |47t 1481 ] given by 



tb 



\r\e 



M 



ir to t/2 



(31) 



for a resonance on top of a background tb ■ For the clastic 
nN residues and their angles as quoted in the PDG, this 
leads to 



\r-KN\ 



arctan 



(32) 



with the arctan function sensitive to the quadrants of the 
complex plane. The phases of the normalized residues 
into inelastic channels are defined accordingly. 

The extracted pole positions, residues, and transition 
branching ratios (normalized residues) are shown in Ta- 
bles IIIII to I VII and Fig. [3H] shows the pole positions in 
the complex plane. Other quantities, usually quoted in 
the PDG, can be trivially calculated from these values, 
in particular the ttN- and other branching ratios, 



r, 



JV 



2\r, 



N 



r, 



r, 



' ^ l/2 ^ l/2 , 

r+ n t 



2|r T 



(33) 



The pN and ttA states can couple to a resonance with 
given J p in more than one way; these possibilities can 
be organized in different channels. The quantum num- 
bers for a value quoted in Tables IIV1 fVl and I VII can be 
obtained by the entry of Table IIXI corresponding to the 
channel number and given J p . 

Quoting transition branching ratios for the 7rA chan- 
nels is physically meaningful because the 7rA state can be 
approximately seen as a pair of stable particles, for ener- 
gies well above the irA threshold. However, as this is not 
the case for the crN and pN channels, we abstain from 
quoting branching ratios for these latter two channels. 
Instead, only the plain residues are shown in the Tables. 
Strictly speaking, the residue of a resonance decaying 
into an on-shell ttttN state is a function of the invariant 
mass of the unstable particle, Mj, (or, equivalently, the 
three- momentum of the stable spectator particle), and 
not just one complex number. For the 7rA state, it is 
a good approximation to fix Mi — Re zq ~ 1232 MeV 
(that kinematics is chosen for the values quoted in the 
tables). 

Note that the labeling of the resonances in Tables IIIII 
to IVII follows the PDG notation event hough the real 
parts of the pole position are different; whether the poles 
determined here can be identified with PDG resonances 
is discussed in the following sections. 

If one defines bran ching ratios using Breit-Wigner pa- 
rameterizations |226| their sum equals, by definition, 1, 
i.e. = I^ot- The right-hand side of this equation 

can be determined independently (r to t = — 2 i Imzo) and 
be used as a test of the present results for the branching 
ratios. Indeed, below the irirN threshold, the equality 
holds at the 1 % level. However, it should be noted that 
= r to t never holds exactly when defining branch- 
ing ratios via residues as done here, even in a manifestly 
unitary coupled-channel model with only stable interme- 
diate states. This is simply because the amplitude has 
non-analytic branch points, required by unitarity, and 
this information is not contained in the residues. 



B. Discussion of specific resonances 

The spectrum of resonances determined in the 
present approach is related to the states found by the 
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TABLE III: Properties of the / = 1/2 resonances: Pole positions zo (r to t defined as -2Imzo), elastic nN residues 
(|fVjv|, 0-xN-iivN), and the normalized residues (yTVivTWrtot, #7riv->/j) of the inelastic reactions nN — > p. For the second 
Roper pole on the hidden nA sheet, shown in brackets, see Sec. IIVD II (*): not identified with PDG name; (a): dynamically 
generated. 
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TABLE IV: Residues of the 7=1/2 resonances for the effective 7r7r# channels: Normalized residues (\/T^Jjl\A/Ftot, GttN^tta) 
for the 7rA channels and coupling constants g for the pN and aN channels according to Eq. (|26|) . For the quantum numbers 
of the ttA channels (6), (7), the pN channels (2), (3), (4), and the aN channel (8), see Table [IX] (*): not identified with PDG 
name; (a): dynamically generated. 
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TABLE V: Properties of the / = 3/2 resonances: Pole positions zo (Ytot denned as -2Im«o), elastic nN residues (|7Vjv|, OnN-tnN), 
and the normalized residues (^/TVivT^/Ttot, O^n^^) of the inelastic reactions nN — > KYI and nN — > nA. For the quantum 
numbers of the irA channels (6), (7), see Table ILXl (a): dynamically generated. 





Pole position 




ttN Residue 




KY channel 




nA, channel (6) 




7rA, channel (7) 




Re zo 




-21m 


zo 






a .. 




pl/2 
1 tvN 


pl/2 
l K E 


OwN- 




pl/2 
1 tvN 


pl/2 
1 irA 


a 




pl/2 
1 tvN 


pl/2 
1 ttA 








|' 7TJV | 




Ttot 


Ttot 




Ttot 




[MeV] 




[MeV] 


[MeV] 


[deg] 




[%] 




[deg] 




[%] 




[deg] 




[%] 




[deg] 


fit^ 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A B 


A(1620) 1/2" 


1600 


1603 


71 


69 


17 


15 


-106 


-107 


21 


22 


-107 


-107 










57 


52 


102 103 




1787 


1848 


575 


379 


54 


19 


-140 


-140 


4.6 


1.9 


-144 


1 1 A 
-114 


39 


22 


71 


80 








A(1232) 3/2+ 


1220 


1218 


86 


96 


44 


50 


-35 


-38 
























A(1600) 3/2+ a) 


1553 


1623 


352 


284 


20 


27 


-158 


-124 


11 


13 


-7 


41 


28 


38 


28 


74 


1 


4 


15 -107 


[A(1920) 3/2+ 


1727 


1813 


866 


892 


31 


19 


169 


-56 


15 


14 


-17 


58 


6 


4 


148 


-98 


1 


<1 


-96 ] 


A(1700) 3/2~ 


1675 


1701 


303 


190 


24 


17 


-9 


11 


1.5 


1.8 


-150 


-112 


5 


4 


167 - 


148 


39 


40 


149 -179 


A(1930) 5/2" 


1776 


1806 


646 


580 


18 


22 


-159 


-103 


3.1 


2.9 


-3 


29 


11 


18 


26 - 


114 


<1 


1 


49 


A(1905) 5/2+ 


1770 


1776 


258 


143 


17 


9 


-59 


-41 


0.5 


0.1 


-142 


-99 


4 


3 


130 - 


179 


34 


29 


105 120 


A(1950) 7/2+ 


1884 


1890 


234 


234 


58 


59 


-25 


-20 


4.0 


3.8 


-77 


-71 


55 


52 


139 


147 


2 


4 


-84 -51 


A(2200) 7/2- 


2148 


2118 


476 


393 


17 


14 


-51 


-41 


0.6 


0.1 


-97 


-75 


2 


8 


-144 


95 


24 


21 


112 121 


A(2400) 9/2- 


1969 


1943 


577 


558 


25 


19 


-80 


-114 


1.3 


0.9 


40 


5 


24 


17 


-98 - 


130 


3 


2 


90 24 



TABLE VI: Residues of the I = 3/2 resonances for the effective pN channels according to Eq. (|26|l . For the quantum numbers 
of the channels (2), (3), (4), see Table ILXl (a): dynamically generated. 



fit^ 


pN, channel (2) 
.g pJ v[10- 3 MeV- 1/2 ] 
A B 


pN, channel (3) 
g pN [10- 3 MeV- 1/2 ] 
A B 


pN, channel (4) 
ffp iv[10- 3 MeV- 1/2 ] 
A B 


A(1620) 1/2" 


4.0-H3.5 3.5+i3.0 






0.9+i0.8 


0.9+i0.9 


A(1910) 1/2+ 


0.8-i2.0 0.4-10.6 


-1.7+i0.9 


-2.2+i3.0 






A(1232) 3/2+ 












A(1600) 3/2+ a) 


0.3-i0.4 1.7-12.0 


-1.9-i3.7 


-1.9-i3.4 


l.l+il.2 


-0.4-il.O 


[A(1920) 3/2+ 


0.8+U.4 2.1+10.9 


4.2+i 6.4 


3.4+i6.5 


0.7+H.2 


0.6+i0.8] 


A(1700) 3/2- 


7.0+i5.7 1.4+i 1.0 


6.3+i6.0 


0.5+H.2 


-6.345.8 


-2.545.7 


A(1930) 5/2 - 


-0. 5-10.2 -0.2+10.4 


1.7+i 1.6 


5.9+i 2.7 


0.8+i0.0 


-0.3+i0.9 


A(1905) 5/2+ 


4.8+i 0.6 3.0-10.8 


4.6+i 1.3 


2.341.2 


-4.0+i 6.3 


-3.5-i0.4 


A(1950) 7/2+ 


0.2-i0.4 -0.1-i0.5 


1.9+i0.3 


-2.8-i 2.6 


0.240.3 


0.3-i0.2 


A(2200) 7/2- 


1.3-11.5 1.241.3 


1.7-i0.2 


1.2+iO.l 


1.4+i5.4 


-0.2+i5.1 


A(2400) 9/2" 


-1.3+i 0.6 -1.0+i0.4 


1.0+i0.1 


0.740.3 


-1.3+i0.5 


-l+i0.4 
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GWU/SAID analysis in Ref. J54J because we fit to the 
corresponding ttN — > ttN partial waves. In that analy- 
sis of clastic and charge- exchange pion-nucleon scatter- 
ing, relatively few resonances are found and the question 
is whether with this set of states, in combination with 
the background provided by the current approach, the 
ttN — > rjN and irN — > KY reactions can be described 
simultaneously with the irN — > irN reaction. It is clear 
that due to our choice of fitting the irN partial waves of 
the GWU/SAID analysis, the most prominent resonances 
of that analysis would be automatically needed here, and 
this is indeed what we find: those are the excited nucleon 
states 



#(1535) 1/2- , 
#(1720) 3/2+, 
7V(1680) 5/2+, 
7V(2220) 9/2+. 



#(1650) 1/2' 
#(1520) 3/2' 
#(2190) 7/2" 



#(1440) 1/2+, 
#(1675) 5/2", 
#(2250) 9/2", 



Note that in Table IIII1 in brackets we have also quoted 
the second pole of the Roper resonance on the hidden it A 
sheet (cf. discussion in Sec UVDljl . 

Furthermore, the following excited A states are found 
both in the GWU/SAID analysis and in the present 
study: 



A(1620) 1/2" 
A(1600) 3/2 H 
A(1905) 5/2 H 



A(1910) 1/2+, 
A(1700) 3/2~ , 
A(1950) 7/2+, 



A(1232) 3/2+, 
A(1930) 5/2~, 
A(2400) 9/2~ . 



The GWU/SAID analysis lists additional poles in the 
complex plane that could not be found in the present 
study, 

#(1860) 5/2+, #(2200) 11/2+ A(2529) 11/2+ . 

But the present analysis includes resonances only up to 
J = 9/2 and thus does not consider the J = 11/2 states, 
and the A(2529) 11/2+ lies outside the considered en- 
ergies anyway. For the second J p = 5/2 + state, see 
Sec. HvUl 

Apart from the states found in the GWU/SAID anal- 
ysis, we include a few more genuine s-channel states that 
are needed for a quantitative description mostly of the 
inelastic reactions irN rjN . KA, and KY,. Those are 



#(1710) 1/2 H 
A(2200) 7/2" 



#(1990) 7/2+, A(1920) 3/2 H 



As Tables UtTI and IV1 show, these states are characterized 
by a weak coupling to the nN channel. 

When discussing pole properties in the following, we al- 
ways refer to the values shown in Tables IIIII to IVII Note 
that branching ratios, when mentioned, can be trivially 
obtained from the normalized residues / transition branch- 
ing ratios in the tables and Eq. (f3"3"| . In the discussion 
of individual resonances, we will highlight particular as- 
pects but we cannot exhaustively compa re ev ery decay 
channel to the corresponding PDG value |226j . 



Sn: The #(1535) 1/2 is quite narrow in both the 
present approach and the GWU/SAID analysis [54| . 
while the #(1650) 1/2" is narrow in the GWU/SAID 
analysis and much wider in the present fits A and B 
(cf. also Fig. ST] in Sec. IIV E[) : still, those fits describe 
the SAID solution quite well on the phy sical axis. In 
the chiral unitary approach of Ref. [45j |. the 7V(1535) 
l/2~ resonance is also narrow [21535 = (l-547+o'o2i — 
i0.046t° n °mi) GcV] while the width of the #(1650) 1/2" 



(1.5971"-- - i0.045™) GeV] 



-0.017 

resonance [zi65o 
lies in between the values of this study and the one 
of the GWU/SAID group. The different values of the 
widths illustrate the difficulties to extract pole positions 
in this partial wave, as will be discussed in Sec. IIV El 
The branching ratio of the #(1535) l/2~ resonance into 
the rjN channel is large (57-59%), while the one of the 
#(1650) l/2~ is much smaller (a roun d 4%), in line with 
what is found in other analyses (226| . The 7r# — > r]N 
residues of the two Sll resonance have a non-zero rela- 
tive phase that indicate a non-trivial interplay, and the 
ttN rjN S-wa,vc cross section shown in Fig. Q3] ex- 
hibits a dip slightly beyond the #(1650) l/2~ resonance 
position; see also Ref. [63( for a discussion. Both Sn res- 
onances show a moderate but non-vanishing branching 
into the 7rA channel. 

The #(1535) l/2~ resonance has a moderate coupling 
to the KA channel and a larger one to the KY chan- 
nel (cf. Table UnJ). The 7V(1650) l/2~ resonance couples 
strongly to all KY channels. Indeed, the partial cross 
sections of the corresponding reactions (Figs. [TH [2TJ l25l) 
show a strong 5- wave contribution close to the thresholds 
and also beyond. It should be stressed however, that 
strong S'-wave components in the KY amplitudes can- 
not be identified with sub-threshold resonances alone, as 
sometimes done in the isobar picture. Instead, one has an 
amplitude with a non-trivial energy dependence in which 
sub-threshold resonances play only a moder ate r ole, as 
demonstrated for pion-induced <p production [l76j . 

P11: For the discussion of the resonance dynamics 
in the Pn partial wave, see Sec. IIVDI We find in the 
current analysis, apart from the (genuine) nucleon and 
the (dynamically generated) two Roper poles, two more 
states: a genuine #(1710) l/2 + state mostly required by 
the inelastic data, and a dynamically generated, broad 
N(1750) l/2 + state, emerging mostly from the cx# and 
7rA dynamics. It should be stressed , however, that the 
#(1710)1/2+ is still around 100 MeV wide and thus can- 
not be identified with a potentially exotic narrow state. 

S31: The width of the A(1620) 1/2" is relatively small 
for both fits A and B which can be attributed to slight de- 
viations of our description from the GWU/SAID partial 
wave, as can be seen in Fig. [7] The resonance shows large 
inelasticities into the ttA channel (normalized residue of 
52-57%, compared to the PDG values of 38± 9%). 

P31: The A(1910) 1/2+ resonance is wide, couples 
only weakly to ttN, and is situated on top of a large, 
energy-dependent background. Due to these intrinsic 
problems, the extracted resonance parameters show large 
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differences between fits A and B. Within these large 
boundaries, the properties agree with those of the PDG, 
in particular with respect to the large inelasticity into 
ttA. 

P13: The 1=1/2 P-wave plays an important role in 
KY production, as, e.g., the A A differential cross sec- 
tion shows. However, there is an interplay between the 
JV(1720) 3/2+ and the JV(1710) 1/2+ resonance and indi- 
vidual contributions are more difficult to pin down. For 
the transition branching ratios [i.e. moduli of the nor- 
malized residues of the A(1720) 3/2+] we obtain the fol- 
lowing numbers (in brackets, the corresponding numbers 
of Anisovich et al [H3, [Hg): rjN: 1.2-3.1% [3 ± 2%], 
A' A: 2.7-3.2% [6 ±4%], KY: 1.8-2%, ttA, P-wave: 22- 
23% [29 ± 8%], and ttA, P-wave: 2-6% [3 ± 3%]. 

Di 3 : The ttN residue of the A(1520) 3/2" state is well 
fixed from ttN scattering where it appears as a sharp and 
prominent resonance. The coupling to rjN is very small, 
but well fixed due to the well-known S'Pz-intcrference 
pattern resulting in the it-shape of the irN — > rjN dif- 
ferential cross section (cf. Fig. [T5")) . We obtain an rjN 
branching ratio of 0.09 — 0.11% while Arndt et al. o btain 
0.08 - 0.12% [H| and Tiator et al. 0.08 ± 0.01% [2271 ]. 
While all these values are comparable, the y diff er from 
the result by Penner et al. of 0.23 ± 0.04% [136j . 

The A(1520) 3/2~ resonance is also characterized by 
a large coupling to the pN channel in S-wave (obviously 
to the low-energy tail of the p meson) , as reflected by the 
very large residue of \g P N\ > 20 in Table HVl There is no 
other state in the -D13 partial wave, although our non- 
resonant part T NP exhibits a resonant structure around 
z ~ 1.7 GcV. This is due to a pole behind the complex 
pN branch point, dynamically generated from the S'-wave 
pN interaction. The dynamics and meaning of this pol e 
are simi lar to the one discussed in Ref. [65|. See Refs. [37l - 
l39l l228j | for further discussion of the dynamical nature of 
this and other states. 

P 33 : Apart from the A(1232) 3/2+, we find a wide 
state that mig ht b e identified with the A(1600) 3/2+ 
from the PDG [226]. As the coupling to the nN channel 
is rather small, compared to the large width, its prop- 
erties are difficult to dete rmin e, and also the pole prop- 
erties cited in the PDG [2261 ] show large uncertainties. 
In our approach, it is dynamically generated from the 
7rA P-wave interaction, showing a significantly larger 
P- wave 7rA residue than in the study by Anisovich et 
al. ^M^M (28-38% vs. 14 ±10%). Note that this state 
was also found in our previous analysis [66|. The second 
genuine resonance state included in the P33 partial wave, 
the A(1920) 3/2+, is characterized by a very high bare 
mass (cf . Tables I VIII IVIII|) that undergoes a large renor- 
malization effect so that the pole moves very far into the 
complex plane. The pole shows a large branching ratio 
into AS and indeed the 7 = 3/2 KY partial cross section 
is large over a wide energy range as Fig. 1261 shows. This 
strength is obviously required by data, although, due to 
the large imaginary part of the pole position, one can no 
longer speak of a resonance signal. 



D33: The A(1700) 3/2~ resonance is relatively wide, 
couples weakly to the nN channel, and is situated on 
a relatively large, energy-dependent background. Out- 
values for the pole position and nN residues agree with 
the PDG [226j within the large intrinsic uncertainties. 

D15: The A(1675) 5/2~ shows a sizeable normalized 
nN ttA residue of 40-50% in line with the PDG val- 
ues [2261 ] (46-50%) but only a small coupling to the KA 
channel. For the reaction nN — > rjN, we find that this 
resonance plays some role as can be seen in Fig. 1141 al- 
though its mass is in a region where the experimental rjN 
data are of particularly poor quality (cf. Sec. IIII D)) . 

F15: The A(1680) 5/2+ resonance shows only small 
branchings into the rjN and KY channels. The role of a 
second P15 state is discussed in Sec lIVCl 

D35: The real part of the pole position of the A(1930) 
5/2~ is in the GWU/SAID partial- wave analysis [Hj 
around 200 MeV higher than in the present study. The 
width is larger in the present fits A and B. In any case, 
the overall size of the partial wave is small and our fit fol- 
lows the trend of the GWU/SAID solution (cf. Figs. 
I47[) : as a consequence, it seems difficult to unambiguously 
determine the pole position in this partial wave, given the 
very large inelasticities into the nnN channels. As can 
be seen in Fig. 1261 the D35 partial wave plays a minor, 
but not insignificant role, in the KY, channel. This is 
reflected in the branching ratio of the A(1930) 5/2~ into 
this channel, which is not large but well determined. 

F35: The A(1905)p35 is a relatively inelastic resonance 
with large branching into nA (normalized nN — > 7rA P- 
wave residue of 29-34% h ere, compared to 25 ± 6% of 
Anisovich et al. [l47l . l226| ). In our result the resonance 
plays no role in KY production. 

F37: The pole position of the prominent A(1950) 7/2+ 
resonance is well determined as Fig. I4T1 shows. This res- 
onance is also very important in the ttN — > KY reaction: 
Our value for the normalized residue and its phase (3.8- 
4% and —77° t o ~71°) is comparable to the value of 
Anisovich et al. [l47|,[226[ (5 ± 1% and -65° ± 25°). 

F17, G37: The resonances included for these partial 
waves mostly serve to improve the description of the nN 
elastic scattering. As Figs. [8] and [48] show, the signal 
of the G37 resonance in irN scattering is weak in the 
sense that the resonance is so wide that the fitted energy 
region covers only part of the resonance shape. The G37 
partial wave is small in KY production. The P17 partial 
wave is relevant for AA production at higher energies 
(cf. Fig. [21]). No resonances are quoted in these partial 
waves in the GWU/SAID analysis [54[ and we cannot 
claim much evidence either. 

G17, G19, H19: The resonances in these partial waves 
are wide, but relatively prominent in elastic nN scatter- 
ing. Our fits to the GWU/SAID analysis [54| are satis- 
factory and the extracted pole positions and n N residues 
are in the vicinity of the GWU/SAID values. The agree- 
ment is better for the G17 resonance and slightly worse 
for the real part of the G19 resonance and the (anyway 
large) imaginary part of the if 19 resonance. None of these 
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FIG. 37: Detail of the J15 partial wave above the prominent 
JV(1680) 5/2 + resonance. Solid (black) lines: energy depen- 
dent GWU/SAID analysis [54J; dashed (red) lines: descrip- 
tion obtained by including a second bare s-channel state in 
the -F15 partial wave and refitting the entire data base (fit A'). 
The faint dash-dotted lines show the original fit A (same as 
in Fig.EJ. 

resonances shows a large branching into the rjN and KY 
channels. 

G39: The G39 resonance is considerably wider in the 
GWU/SAID partial-wave analysis than here (900 MeV 
vs. 550 MeV), but for poles so far in the complex plane, 
one expects uncertainties to become very large. The nor- 
malized residue into KY, is small (0.9-1.3%) but seems to 
be significant, as suggested by the relatively small uncer- 
tainty for the partial cross section in Fig. l26l 



C. Fit A': a second F15 state 

Coming back to the -F15 partial wave, only one gen- 
uine resonance is included in fits A and B, while two 
resonances are needed in the GWU/SAID analysis. The 
second F15 resonance state, iV(1860) 5/2 + in the PDG, 
formerly known as V(2000) 5/2+ , is not included in the 
present analysis; from Figs. |5] and 25] it is obvious that 
a more quantitative description of the F15 may require 
this state. As a test, we show in Fig. [37] the result when 
including a second bare s-channel state in the -F15 par- 
tial wave and performing a refit (fit A') of the entire 
data base, starting from fit A and varying all bare reso- 
nance masses and couplings. The description of the -F15 
partial wave around z ~ 2 GeV improves; in fact we 
achieve an almost perfect match of the GWU/SAID so- 
lution. However, we see no substantial improvement for 
the other reactions with rjN and KY final states. The 
largest changes are visible in the KA polarization data 
at backward angles, but the data carry large errors (cf. 
Fig. [TTJ. Thus, the inelastic reactions cannot contribute 
much to increase the evidence for this state, and more 
precise polarization data are called for. 

Even though for fit A' we have varied all resonance 
parameters, the properties of other states change only 
marginally. Also, the properties of the strong V(1680) 
5/2 + resonance barely change comparing fits A and A' 
[20 (fit A) = (1666 - m.5) MeV vs. z (fit A') = (1670 - 



i58) MeV]. In fit A', we obtain \r nN \ = 40 MeV compared 
to \r nN \ = 42 MeV in the SAID solution 0. The pole 
position of the second F15 resonance AT(1860) 5/2 + is 
at 2 (fit A') = (1730 - i66.5) MeV, while the residue 
amounts to |jViv| = 7 MeV and O-nN-urN = 180°, i.e. is 
much smaller than the one of the V(1680) 5/2 + . This is 
in line with what is naively expected from the shape of 
the amplitude, i.e the iV(1860) 5/2 + resonance appearing 
as a small bump to the right of the prominent AT(1680) 
5/2+. In Ref. in contrast, the V(1860) 5/2+ has a 
residue of Ir^jv | = 60 MeV. This value, even larger than 
the one of the JV(1680) 5/2 + , might have to do with a 
close-by zero and the spec ial dynamics of this resonance 
in the SAID analysis |229j . 

D. The Pn partial wave in ttN — > nN 

1. Fits A, B, and the Roper resonance 

The P\i partial wave and its analytic structure require 
a detailed analysis. As mentioned in Sec. IIII Bl we use 
different starting conditions for fit A and B. In particular, 
we demand the nucleon renormalization to be small in 
fit B. As a result we observe, in the course of fitting, 
differences in the dynamics for the inelasticities in the Pn 
partial wave, which arc closely linked to the properties 
of the aN channel. In the following, we illustrate the 
interplay between nucleon renormalization and changes 
in the aN channel. 

The aN channel has much influence in the Pn partial 
wave due to the onset of large inelasticities at relatively 
low energies and the fact that in Pn the two pions and 
the nucleon of the aN channel are all in relative S- waves, 
i.e. there is no centrifugal barrier to suppress low-energy 
contributions. 

By far the largest differences in the parameters found 
for fits A and B are those in the g a NN and g aaa coupling 
constants, cf. Table IXlTTl The implications of a small or 
large g a NN, or generally speaking of the two scenarios A 
and B, is illustrated in Fig. 1381 where we show the non- 
pole part of the r matrix for both cases [r NP and r p are 
related to T NP and T p via Eq. {2"7|]. 

For fit B (large g a NN), the pole in r NP is obviously 
already in the vicinity of the full amplitude r. The pole 
part, r p , consists of the nucleon plus a second pole re- 
sponsible for the V(1710) l/2 + that we can omit from 
the discussion for the moment. When adding the nucleon 
pole, it repels the pole in r NP due to the mechanism of 
resonance repulsion [&4| and shifts it to its final position 
in the full amplitude r = t p + r NP . As this shift is 
rather small for fit B, so is the effect due to renormaliza- 
tion of the nucleon pole. Indeed, the bare nucleon mass 
of m b = 974 MeV (cf. Tabic |VTTT|) is close to the phys- 
ical one. Also, the renormalization of the bare irNN 
coupling (/^, = 1.020) is small when compared to the 
dressed, physical value of f^NN = 0.989. 

For fit A (small g CT jvAf), the scenario is different: here, 
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FIG. 38: Imaginary part of the Pn irN amplitude. Besides 
the full amplitudes r for fit A (solid black line) and fit B 
(dashed black lines), also the corresponding non-pole parts 
r NP according to Eqs. ((6] I27|) are shown. Note that this de- 
composition has no physical meaning and serves only for illus- 
tration of the renormalization mechanism due to the nucleon 
pole term. 
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FIG. 39: Analytic structure of the present fit result (fit B) 
in the Pn partial wave, for complex scattering energies z. 
The circular cut is labeled "CC", the left-hand cut "LHC". 
See text for further explanations. Note that in fit A the bare 
mass of the nucleon is above the irN threshold. 



In summary, while the internal parameterization of the 
Pn partial wave is quite different, the full, physical am- 
plitude is not, and in both scenarios, it is possible to 
quantitatively describe the GWU/SAID amplitude. 



the pole in r is farther away from the Roper pole po- 
sition as Fig. [38] shows. Due to the discussed mechanism 
of resonance repulsion, the nucleon renormalization effect 
is larger and indeed we observe a bare nucleon mass of 
m b = 1145 MeV (cf. Table EH]). The large renormaliza- 
tion of the bare mass in fit A comes along with a large 
renormalization of the bare coupling that is = 0.646 
compared to the physical, dressed value. The size of the 
renormalization of fit A is comparable to the one of pre- 
vious studies [62], HII- For both scenarios A and B, we 
obtain a slightly negative phase shift close to the it N 
threshold (cf. Fig. |4"5|) which is a remnant of the nucleon 
bound state and its interplay with the Roper resonance. 

While now the renormalization of the nucleon is under- 
stood for fit A and B, the question remains what further 
differences there are in the parameterization of the T NP 
amplitude. Tracing back the origin of the pole in T NP , 
it turns out that for both fits (and the previous ones of 
Refs. [59 63]) a pole in the Pn partial wave is generated 
from the ^-channel correlated two-pion exchange with p 
quantum numbers. For fit A (small g^NN), adding all 
other t- and w-channel exchanges, one ends up with the 
pole structure visible in Fig. [35] For fit B (large g a NN), 
the exchange contributions to and within the aN channel 
— see Figs. [5] and [3] — provide so much attraction, that the 
pole generated from the p-exchange is driven far below 
the nucleon pole position and becomes a spurious state 
on one of the many Riemann sheets there (cf . Fig. l39l) . 
The physical amplitude below the nucleon pole and be- 
tween nucleon pole and ttN threshold is pole-free. The 
large aN attraction provides at the same time enough 
strength for the formation of the pole visible in T NP of 
Fig.[35J 



2. Analytic structure 

The analytic structure of the current fits A and B is dis- 
played in Fig. [35] Shown are the branch points (thresh- 
olds) of the ttN, r)N, KA, and KYj channels on the real, 
physical z-axis, as well as the ttttA^ branch points on the 
real axis and in the complex plane from the 7rA, aN, 
and pN channels. To search for poles on the most rele- 
vant Riemann sheets, the cuts are chosen to be oriented 
from the threshold openings into the negative imaginary 
z-direction as indicated in the figure (see also Sec. IIV Al) . 
In addition, the figure shows the sub-threshold short nu- 
cleon cut (SNC), circular cut (CC) and left-hand cut 
(LHC). The structure of cuts and branch points is dis- 
cussed in detail in Ref. [65| . 

The (green) stars indicate the positions of the physical 
poles. There are two Roper poles on different 7rA sheets, 
labeled as N(1440) and N'(1440) with the positions given 
in Table IIIII This two-pole structure is found in several 
analyses, such as GWU/SAID H and EBAC (z^, and is 
due to the fact that poles can repeat on different sheets. 
The closeness of the complex ttA branch point to the 
resonance poles leads to a non-standard shape that is 
not of the Breit-Wigner type. 

Below the irN threshold, the cross indicates the bare 
nucleon pole and the (green) star the physical position 
of the nucleon. The same labeling is used for the second 
genuine state [the iV(1710) l/2 + ] with the pole positions 
given in Table IIIII The properties of this state are dis- 
cussed in the following section. In the irN amplitude 
shown in Figs. 151 and P4"5l it is visible as a weak resonance 
above the Roper state. 

When comparing the pole properties of the Roper res- 
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onance for fits A and B according to Tables Mil and IIV1 
it turns out that for fit B the position is considerably 
farther in the complex plane and the residue for decay 
into the aN channel is smaller than for fit A. In turn, 
the residue for decay into the irA channel is larger for fit 
B. This demonstrates that the net inelasticity, that must 
be approximately the same for fit A and B, is distributed 
differently among the effective ttttN channels for the two 
fits. 

Finally, we find one additional, dynamically generated 
pole in the complex plane at around z = 1.75 GeV and far 
in the complex plane, called N(1750) 1/2+ in Tables Hill 
and|lVl This state appears in both fits A and B, but with 
a very large aN residue in fit B and a small one in fit A. 
We expect that a dedicated study of the reaction ttN 
ttttN can shed light on the question which of the two 
scenarios A or B provides the more realistic description. 

Here, we just note that the analytic structure of Fig. 131)1 
is different from the one found by the EBAC group in 
Ref. [78[ : in that parameterization, there is only one gen- 
uine resonance state with a bare mass at z = 1736 MeV, 
that, through the hadronic dressing, is responsible for the 
two Roper poles plus a pole at z = (1820 - i248) MeV. 
The nucleon pole is not included in Ref. [78^ . however, in 
Ref. [22l ! the EBAC group investigated the consequences 
of a genuine nucleon pole. 

The question arises if there are ways to distinguish the 
reaction dynamics found here from the one by the EBAC 
group in Ref. [78]. The principal problem is that the 
only observable quantity is the fully dressed amplitude, 
while quantities such as bare masses and couplings de- 
pend on the renormalization scheme and the space of in- 
cluded channels, as demonstrated in Sec. lIVFl Therefore, 
while certain scenarios might be ruled out in a dedicated 
study of 7T7rA final states, there is, in principle, no way 
to distinguish between internal parameterizations of the 
amplitude in terms of bare quantities, and, in particular, 
no meaningful hadronic "undressing" of resonances. 



3. Properties of the jV(1710) 1/2^ 

In Fig. l40l the region around z ~ 1.7 GeV is displayed 
in more detail. The figure shows, apart from fits A and 
B, the energy-dependent Pn partial wave from Ref. [54[ 
that we fit (green dash-dotted lines). Most intriguingly, 
our fit B, and to some extent also A, match the single- 
energy solution from Ref. [54[ better than the energy- 
dependent solution (green curve) we fit — the latter is 
smooth around z = 1650 MeV. 

In the present analysis, the inclusion of this second 
resonance became necessary to improve the ir~p — > K°A 
differential cross section and polarization, and was not 
required at all for the fit to the energy- dependent irN 
SAID solution. Note that in the fit procedure, the reso- 
nance position was not forced in any way to its current 
position, but left to float freely. 

Thus, the structure seen in the single-energy solution of 
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FIG. 40: Detail of the partial- wave amplitude (top) and phase 
shift/inelasticity (bottom) of the elastic nN — s> nN P\\ wave. 
The current results, fit A (solid red lines) and fit B (dashed 
blue lines) , were fitted to the energy-dependent GWU/SAID 
solution [541] (dash-dotted green lines). The single-energy so- 
lution [5J| (points) is not included in the fitting procedure. 

the analysis of ttN —> ttN [54| appears here for an entirely 
different reason, namely through coupled-channel effects 
from the KA data and, to some weaker extent, the r/N 
data. We see that as evidence for a resonance at this 
energy. 

In Ref. [8!|, the role of the pN branch point, or com- 
plex threshold opening, at z ~ (1-7— i 0.07) GeV, i.e. 
not so far away from the discussed resonance, was scru- 
tinized. It is shown that if an approach does not contain 
such branch points, false resonance signals may appear. 
However, the current approach does contain the required 
pN, aN, and ttA branch points, and still we observe a 
new resonance around z ~ 1.65 GeV. We note that the 
resonance has a considerable branching fraction into r/N 
and KA; cf. Table [IH] See also Figs. Q3 and [21] for the 
partial-wave content of the ir~p — > T]n and n~p — > K°A 
reactions, respectively. For the r/n final state, the Pn 
resonance is partly responsible for the second maximum 
of the total cross section, although it should be stressed 
that the data situation in that energy range is very poor; 
see the discussion in Sec. IIII Dl For the A' A final state, 
Fig. [2~T1 shows that in our fit the strong P-wave contribu- 
tion at z ~ 1.7 GeV comes from the Pn wave instead of 
the Pi 3 wave that could also play a role here. This is ac- 
tually demanded mostly by the polarization data in that 
energy range. Here, the partial- wave content is fixed bet- 
ter than in the rjn case because the data are much better. 

The role of a Pn resonance at z ~ 1.7 GeV in 
the context of different partia l- wav e analyses has been 
discussed extensively in Ref. |152| (see also references 
therein). That resonance appears in several quark mod- 
els fi~7T4l£l but also in three-body hadronic calcula- 
tions [230| . In Ref. [l52j . various fits to the single-energy 
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FIG. 41: Pole positions of the 7V(1650) 1/2" (left) and the 
A[1950) 7/2+ (right) in different analyses [II, H El SI 11 
1 146f | . The results of this study are indicated with the (red) 
square (fit A) and (blue) diamond (fit B). Also, the KA, KT,, 
and pN branch points are shown (green crosses). 



GWU/SAID solution and existing r]N partial waves were 
performed and the need for an V(1710) l/2 + resonance 
was stated. Here, in contrast, we fit directly and si- 
multaneously to the r]N and KA data and the energy- 
dependent GWU/SAID solution, i.e. the present fit is 
unbiased by structures in the single-energy GWU/SAID 
solution ("data" points in Fig. l40]) . Also, instead of being 
required by the (rather low-quality) rjN data, the need 
for a P\i state around z ~ 1.7 GeV arises from the KA 
data. In any case, it is noteworthy that in two different 
analyses the need for aPn state around z ~ 1.7 GeV is 
found. 



In Ref. [23 ll ] , various ttN partial waves were scanned in 
the search for narrow resonances. Two candidates were 
found; a stronger signal at z ~ 1.68 GeV and a weaker 
one at z ~ 1.73 GeV in the Pn partial wave. How- 
ever, while that first possibility would be in line with the 
current findings, the resulting widths were much smaller 
than the one found here, which is around 120 MeV (cf. 
Table IIIip . We do not make any attempt to identify the 
resonance found here with a narrow state of potentially 
exotic nature; it is just worth noting that the present ap- 
proach does contain the KY, channel that has its thresh- 
old close to the observed structure; threshold effects from 
this channel, are, thus, taken into account. Such thresh- 
old effects can become large in S'-wave and might be even 
responsible for the excess of r\ photoproduction on the 
neutron around these en ergies [232| (for other explana- 
tions, see Refs. (233T - l235| ); we do not observe any struc- 
tures from KT, in the Pn partial wave in the current 
fit result. Analyzing KY photoproduction data in the 
future might shed more light on the issue. 
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FIG. 42: Present result (fit A, solid lines) together with the 
Laurent expansions of Eq. (|34[1 around the different poles in 
the Sn partial wave (dashed lines) and the sum of both terms 
(dashed-dotted lines). 



E. Uncertainties in pole positions and bare 
parameters 

In Fig. [4TJ the present results for the pole positions 
of the iV(1650) 1/2 - and the A(1950) 7/2+ are com- 
pared to other analyses quoted by the PDG j226| . For 
the isolated F-wave resonance, the results of the different 
analyses agree quite well. By contrast, the results for the 
V(1650) 1/2" S'-wave resonance vary much more, be- 
cause there is another close-by state, the iV(1535) 1/2 - . 
To illustrate this effect, in Fig. [42] we show the present 
description of the Su partial wave (fit A, red solid lines) 
together with the quantity 



Z - Z 



(34) 



for both Six resonances (green dashed lines) where tvjv 
is the elastic residue and zq the pole position of the re- 
spective resonance. Obviously, the tail of one resonance 
provides a strongly energy-dependent background for the 
other resonance. Indeed, by summing the contribution 
from both resonances, one obtains the dash-dotted curves 
that resemble the full amplitude quite well. 

However, the mutual influence of close-by resonances 
makes their extraction challenging. As Table UTT1 shows, 
in this study we obtain a small width for the iV(1535) 
l/2~ in combination with a rather large width for the 
vV(1650) 1/2". Note that also the GWU/SAID group 
found a relatively narrow width for the iV(1535) 1/2 - 
[54| in 2006, although the resonance width determined in 
2004 [53l is larger by 40%, which again illustrates the dif- 
ficulties in resonance extraction for the S\\ partial wave. 

Despite the fact that the extraction of the S\\ res- 
onances is, for the above reasons, certainly one of the 
more complicated cases, the real parts of the pole posi- 
tions differ only by 2 MeV in our approach (cf. Tab. IIII[) 
whereas the bare mass are separated by around 50 MeV 
(cf. Tab. EU] and Gnu} . 

Increased uncertainties for bare masses, compared to 
pole positions, can be also found for other resonances 
comparing fits A and B. For example, the real parts of 
the pole positions of the iV(1710) l/2 + resonance, that 
by itself is difficult to pin down, differ only by 18 MeV 
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(cf. Tab. lIII[) while the bare masses are apart by 143 McV 
(cf. Tab. IVlil and IVTIl]) . Similar considerations hold for 
the very prominent A(1950) 7/2 + resonance (differences 
of 6 MeV in Re z vs. 98 MeV in m b ). 

This illustrates the high uncertainties tied to bare pa- 
rameters, cf. also the discussion at the end of Sec. Ill Bl 
In the following section, we will further discuss this issue. 



F. Residues vs. bare couplings 

To point out the large correlations between bare reso- 
nance couplings and bare masses, we consider the case of 
the A(1620) l/2~ resonance, characterized in the present 
fits by very high bare masses. For the test, starting from 
fit A, we reduce the bare mass m b from its actual value 
of 4.8 GeV (cf. Tab. EH]) and then perform a refit of 
the bare resonance couplings, i.e f nN , f pN , f n A, Iky,-, to 
the full data base. For different bare masses, we show in 
the upper panel of Fig. 03] the results for the amplitude. 
Reducing the bare mass, the quality of the refits is in- 
distinguishable down to m b = 2.9 GeV (changes in the 
X 2 are less than 1%). When the bare mass is further re- 
duced, the amplitude shows noticeable deviations. In any 
case, the real part of the pole position is shifted by 2 McV 
at m b — 2.4 GeV compared to the case m b = 4.8 GeV. 
The imaginary part is shifted by —5 McV. This means 
that changes in the bare mass can be almost fully com- 
pensated by changes in the bare couplings. 

In the middle panel of Fig. 23J we show the relative 
changes of the bare couplings as m b is changed (filled 
symbols). Also, for the case m b = 2.4 GeV, the relative 
changes of the residues are displayed (for this mass, the 
X 2 has changed by 2.5%). While bare couplings change 
by up to 60%, the coupling constant g in the ttN channel 
[i.e. the square root of the residue, cf. Eq. ([2"o])] changes 
by 4% (empty square). The coupling g of the 7rA channel 
changes by 10% (empty triangle down), and the one of 
the KH channel by 17% (empty triangle right). The size 
of these changes is expected, given that the irN ampli- 
tude is precisely known while the ttA channel has only 
indirect influence through the resonance width, and the 
KT, channel is still closed at these energies (the resonance 
still contributes to the description of KT, data, and thus 
to the x 2 , through its finite width). 

In the lower panel of Fig. 23] the contribution of the 
different channels to the self energy £ is shown accord- 
ing to Eq. Q, i.e. every term in the sum over the 
channel index fj.. The self energy is evaluated at the 
pole position. As the figure shows, the real parts of the 
self energies are very different for different bare masses, 
fully compensating changes in the bare mass to ensure 
an almost unchanged pole position. For the imaginary 
part, shown to the right, one should note that differ- 
ent channels contribute with different sign. As indi- 
cated with the horizontal black lines, the net half width 
Ttot/2 = —Ihizq = — Im£(z = zq) is rather small and 
does not change significantly for different bare masses 
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FIG. 43: Upper panel: Dependence of the 531 amplitude on 
the bare mass of the A (1620) l/2~ (involves refits of the four 
bare couplings /(,). The case m b — 4.8 GeV corresponds to 
fit A. Middle panel: change of the /;,, normalized to the case 
m b = 4.8 GeV (filled symbols). For m b = 2.4 GeV, also the 
changes of g from the residues [cf. Eq. (|26p ] are shown (empty 
symbols). Lower panel: contributions to the self energy E(z = 
zo) from different channels. 



(by 5 MeV) although the contributions from individual 
channels change dramatically. 

This numerical example again illustrates the problems 
tied to bare parameters. On top of this, there are also 
large correlations between bare parameters and the cut- 
off in the form factors of the resonance vertices that we 
have not investigated for this example. It is, however, 
clear that a change in the cut-off can be compensated 
by one in the bare coupling; thus, there is an additional 
scheme dependence of the bare coupling constants (and, 
of course, also of the bare mass through the self-energy), 
that further prohibits to attach a physical meaning to 
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these model parameters. 



G. Outlook: significance of resonance signals and 
error analysis 

The presented fits A and B are the result of exten- 
sive searches in parameter space, performed under dif- 
ferent starting conditions to ensure that in the vicinity 
there are no better minima. Nonetheless, rather than 
having several well-distinguished local minima of the % 2 , 
one observes many minima that are shallow in the sense 
that there are strong parameter correlations, for example 
between bare resonance masses and coupling constants 
as discussed in Sec. IIVFI In general, if the contribution 
to the imaginary part of the resonance self-energy [cf. 
Eq. ||3J] is small, as in case for the pN channel, a change 
in the bare mass can be absorbed well by a change of 
the bare coupling. Quantifying such correlations is post- 
poned to a later stage of the analysis, but is, in principle, 
feasible within the current approach, as has been shown 
in Ref. @. 

As for the statistical error (in contrast to the system- 
atic error as estimated here by performing two different 
fits), in Ref. (6(| it has also been demonstrated how the 
error of the experimental data translates into uncertain- 
ties of extracted resonance pole positions. It should be 
noted, however, that a true statistical analysis will be 
only possible once one fits to the actual ttN — > ttN data 
and not, as was done here and in most other currently 
active analysis efforts, to the partial- wave solution of the 
GWU/SAID group. We have fitted here to their energy- 
dependent solution as recommended by the GWU/SAID 
group, but it should be also noted that the error bars 
provided in their energy-independent (SES) solution do 
not have statistical meaning, i.e. they cannot be used 
for the purpose of statistical error analysis in a multi- 
channel approach of the present type. This is also one of 
the reasons why we refrain from quoting an overall \ 2 ■ 

In addition to studying the propagation of the statisti- 
cal error from the data to results such as the partial- wave 
decomposition and pole properties, there are systematic 
uncertainties tied to the selection of experimental data 
in the fit and the weight of those data in the fitting pro- 
cedure. As discussed before, the very different quality 
of the data makes it unavoidable that some reactions 
and data sets are weighted differently from others (apart 
from different weighting due to additional systematic un- 
certainties in experiments, as discussed in the context of 
particular reactions in previous sections). In other words, 
the present fits are, unavoidably, the result of a selection 
process. 

Aside from fitting directly to irN data, another pos- 
sibility to improve the extraction of the resonance con- 
tent is to include photon-induced reactions. A big step 
in this direction has been recently achieved by fitting 
pion photoproduction data using the gauge-invariant ex- 
tension of the present framework [TOj . With the in- 



clusion of KY channels in the present study and the 
good description found for pion-induccd 77./V, KA, and 
KY, production, it becomes now possible to analyze also 
the high-precision data on photoproduction of these final 
states. Determining the photocouplings to excited bary- 
onic states and their Q 2 dependence is also important 
for other th eoret i cal appro aches that need analyses such 
as in Rcfs. [l42l 1 2 3 6l-l 2 3 8|j as input, namely chiral uni- 
tary calculations [33ll35l l45l 123914241 1 that are sometimes 



compared to multipolcs, in particular for pion photopro- 
duction. A simultaneous fit to both pion- and photon- 
induced reactions will increase the reliability of resonance 
extraction possible in the current approach due to its cor- 
rect analytic properties. 



V. CONCLUSIONS 

A combined analysis of the reactions irN — > irN, rjN, 
KA, and the three measured KY final states K + Y + , 
K°Y°, and K + Y~ within a dynamical coupled-channel 
framework has been performed. For the inelastic re- 
actions, the world data set from threshold to ener- 
gies around 2.3 GeV has been considered. In the 
Lagrangian-based calculation, the full off-shell solution 
of the Lippmann-Schwinger-type equation provides the 
correct analytic structure allowing for a reliable contin- 
uation into the complex plane to extract resonance pole 
positions and residues up to J p = 9/2 ± . The amplitude 
features also the effective ttttN channels 7rA, pN , and aN 
with branch points in the complex plane and a dispersive 
treatment of a and p t-channel exchanges. 

The Lagrangian-based, field-theoretical framework is 
particularly suited to perform coupled-channel analyses: 
the SU(3) flavor symmetry for the exchange processes al- 
lows one to relate different final states, and the t- and u- 
channel diagrams, whose explicit inclusion is mandated 
by three-body unitarity arguments, connect also differ- 
ent partial waves and the respective backgrounds. As a 
result, for all reactions, a realistic and structured back- 
ground can be provided that helps map out resonance 
signals. 

Systematic uncertainties have been roughly estimated 
by performing two different fits with starting points far 
from each other in the parameter space. Although the 
parameterization of the amplitudes is quite different, pole 
positions and residues change much less. In contrast to 
these pole parameters, bare masses and couplings were 
shown to have no physical meaning. 

All of the well-established four-star resonances within 
the considered energy range were found and their branch- 
ing into the rjN and KY channels determined. For 
resonances less prominently visible in irN scattering, a 
P\i state of around 120 MeV width was found, required 
mostly by the KA data and coupled-channel constraints. 
Indeed, the state might be identified with the ^(1710) 
l/2 + . A further, but rather wide state, has been found in 
this partial wave. Furthermore, from the KY, data there 
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is need for large strength in the P33 partial wave around 
z ~ 1.9 GeV. A A(1600) 3/2+ state found in a previous 
analysis could be confirmed. 

Apart from determining the resonance spectrum, also 
scattering lengths and volumes have been calculated 
which are in fair agreement with approaches using chi- 
ral perturbation theory; the consequences of adding and 
removing a second F15 state has been tested and the im- 
portance to take the data of all three K E final states into 
account has been demonstrated. 

The present results, in combination with the recent ex- 
tension to pion photoproduction [70l ] , constitute a major 
step towards the analysis of high-precision photoproduc- 
tion data of 77 AT, A'A, and ICE data produced, e.g., at 
ELSA, JLab, and MAMI, allowing for a more fundamen- 
tal understanding of the resonance spectrum. 
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Appendix A: Bare resonance parameters 

In Tables IVIII and I Villi the bare resonance parame- 
ters are quoted for reproducibility of the fit result. As 
explained in Sec. Ill Bt bare masses and coupling con- 
stants have no physical meaning, because they depend on 
the scheme (i.e., values of the form factors), the number 
of included channels, and some bare resonance param- 
eters tend to be strongly correlated (cf. Sec. \TV F|) . In 
fact, some bare masses were constrained by hand because 
large values were obtained in the fit, without noticeable 
improvement of the \ 2 compared to lower bare masses. 
Thus, we arc aware that some masses and couplings in 
Tables FTTH and IVIIII arc rather large, but, considering the 
above arguments, we regard those high values as unprob- 
lcmatic. We expect that once ttttN data are included in 
the fit, the bare masses and couplings tied to the ttttN 
channel can be fixed more reliably (their values still being 
scheme-dependent), but at the current state of the ap- 
proach it makes no sense to limit these values too much. 

The coupling constants / of Table IVIII and IVIIII enter 
the vertices that couple the JV*'s and A*'s to the dif- 
ferent channels. The corresponding Lagrangians can be 



found in Table 8 of Ref. [66|. The partial- wave projected 
vertex functions, for vertices up to J = 3/2, are given in 
Appendix B of Ref. [gf| . In this study, we also need the 
vertices of J = 9/2 resonances. For this, we use the pre- 
scription analogous to the J = 5/2 and J = 7/2 vertices, 
given in Eq. (B.3) of Ref. (66|; i.e. the J = 9/2 vertex 
functions are obtained from the J = 3/2 vertex functions 
according to 

(7 a 'V =$(7-V 

(7 a 'V =^(7 a 'V (Al) 

with the functions (7 a ' c )a± from Eq. (B.2) of Ref. 
These 7 a,c contain the partial-wave projected vertex 
functions v, isospin factors Ir (cf. Table 10 of Ref. j66j|). 
and form factors F with a common cut-off of A = 2 GeV, 
except for the nucleon cut-off, which is a free parameter 
in our model and determined to be A^ = 1575 McV in fit 
A and A N = 958 McV in fit B. Equation ((ATI) provides 
the correct dependence on the orbital angular momen- 
tum L for all channels MB = nN, ?]N, KA, KT,, pN, 
and 7rA; see also Table IIXI In principle, it would also 
be possible to derive the vertices for J > 5/2 f rom La- 
grangians, as has been done, e.g., in Ref. [139j for the 
ttN channel. 

Note that there arc different possibilities for the 7rA 
and pN channels to couple to a given J p . For conve- 
nience, we show the complete coupling scheme (applying 
to s-, t-, and u-channcl exchanges) in Table IIXI There, 
the total spin S = \Sn + S p \ is given by the sum of the 
p spin and the nucleon spin, and L is the orbital angular 
momentum. The notation of the quantum numbers is the 
usual one, L>2i-,2J, e.g. Su, etc. Different coupling pos- 
sibilities for the pN and 7rA channels are incorporated 
here as separate channels p. 

Appendix B: Exchange potentials 

In Appendix IB II we list explicit expressions V for all 
t- and u-channel exchanges used in this work and shown 
in Figs. [2] and [3] The corresponding interaction terms of 
the Lagrangians can be found in Table [X] Vertices that 
are not listed in Table [X] are related to the ones in the 
table by SU(3) flavor symmetry. The explicit expressions 
g a , gb for those couplings can be found in Sec. IB 21 Note 
that for the /o exchange the same Lagrangians as for the 
exchange of a stable a were applied. 

In the following, we give further details on the nota- 
tion used in in Appendix IB II The labeling of particles 
and helicities are specified in Fig. [44] The index 1 and 
3 (2 and 4) denote the incoming and outgoing baryon 
(meson). The on-shell energies arc 

Ei = ^P?+mj, u t = yjtf + m? (Bl) 

for the baryon and the meson, respectively, and mi, 
m 3 ( m 2, "14) are the masses of the incoming, outgoing 
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TABLE VII: Bare resonance parameters fit A: masses m b and coupling constants /. Note that the bare couplings differ 
sometimes by orders of magnitude among each other, because the power of the form factors can be different to ensure convergence 
of the integrals. In particular for the high partial waves, those powers can be large. 
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baryons (mesons). In the TOPT framework used in this 
study, the zeroth component of the initial and final mo- 
menta are set to their on-mass-shell values: p® = Ei or 
Pi = Wi. 

The three-momentum of the intermediate particle is 
q. The four-vector q with q° — E q (baryon exchange) 
or q° = ojq (meson exchange) means the exchanged mo- 
mentum in the first time ordering whereas q indicates the 
second time ordering with q° = —E q (baryon exchange) 
or q° = —u! q (meson exchange). For the potentials, we 
use the notation f = 7 M p M . 

Each exchange diagram includes a kinematic normal- 
ization factor 



N = 



1 



1 



(2tt) 3 2^7^ 



(B2) 



Isospin factors (IF) for the exchange processes (we use 
the baryon-first convention) can be found in Tables IXII 
and IXIIl Also, every exchange process V quoted in Ap- 
pendix IB II is multiplied with form factors Fi , F^ that 
depend on the exchanged momentum q, quoted in Ap- 
pendix IB 31 and Table IXIVI The final expression for the 
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FIG. 44: t- and u-channel exchange processes. 



t- and u-exchanges is then given by 



v^:sis:=n Fi f 2 (if)v 



(B3) 
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TABLE VIII: Bare resonance parameters fit B: masses m b and coupling constants /. 
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TABLE IX: Angular momentum structure of the coupled channels in isospin / = 1/2 up to J = 9/2. The I = 3/2 sector is 
similar up to obvious isospin selection rules. 
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To obtain from this the transition potentials V in particle with momentum pi and helicity A,. 
Eq. (|4]) in the JLS basis, the projection techniques de- 
scribed in Refs. Ell] are used. 

With m ex we denote the mass of the exchange particle. Furthermore, in the potentials quoted in Appendix lB 11 

t v {pi,\i) is the polarization vector of a massive spin 1 the Rarita-Schwingcr spinor for spin 3/2 particles is given 
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TABLE X: Interaction terms of the effective Lagrangians. We use p^ v — <9 M p„ — <9„p M . The vector notation refers to isospin 
space. 
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and is the Rarita-Schwinger tensor of spin 3/2 par- 
ticles (el, 



In case of vector mesons coupling to octet baryons, 
Ibbv = 9bbv denotes the tensor coupling, gsBV 
is the vector coupling (see Appendix IB 2[) ; e a $nv is the 
totally antisymmetric Levi-Civita tensor with e a| g^„ = 1 
for odd permutations of a, (3, p, v, and a^ v = i( 7 M 7 I/ — 
7"7")/2. 



U V 

. 2 " 

3 m£ x om ex 



W-f-f-f)] ■ (B5) 



If pv v appears with the argument qs* for the exchange of 
a J p = 3/2+ baryon B* = A(1232), E*(1385), S*(1530) 
in Appendix IB 11 we use the simplified prescription from 
Ref. [60( in order to avoid problems arising from the con- 
struction of the contact graphs in TOPT (60[. In that 
case, the zeroth component of the exchanged momentum 



is q° B 
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E° n = 



z 2 + m? 



2 z 



z 2 — m§ + m\ 
2z 



(B6) 



Note that for the exchange of p and a quantum num- 
bers in the ttN — > ttN transitions, we use correlated 
two-pion exchange. By contrast, for the exchange of 
these quantum numbers involving other channels, no 
pseudo-data exist to fix these transitions (cf. Sec. [TTJ) 
and we resort to the exchange of stable particles with 
m p = 770 McV and m a — 650 McV. For other t-channel 
processes, the exchanged masses are = 783 MeV, 
m fo = 974 MeV, m ao = 983 MeV, m ai = 1260 McV, 
m K * = 892 McV, and m = 1020 McV. The masses of 
exchanged baryons take their PDG values [226l |; for the 
unstable E*(1530), we use m H * = 1533 MeV. 



1. Amplitudes for the exchange diagrams 

Here, we define the exchange pseudo-potentials V. Common factors have been isolated according to Eq. (|B3I) . The 
amplitude for the correlated mr exchange in the p channel (cf. also Fig. 2]) can be found in Eq. (A7) of Ref. [62| . For 
the correlated tttt exchange in the a channel, see Ref. (63|, Eq. (A3). For this transition, we use the same subtraction 
constant of A = 25 MeV/F 2 as in Ref. (63| . as we see no need to change its value. 

Below, a "type" number is introduced for a more compact notation. To obtain the amplitude for a specific diagram 
(cf. Figs. [5] and [3]), one uses Eq. (|B3[) with the form factors F\, i*2 from Appendix IB 31 V from this section according 
to Tabs. IXII and IXIII and isospin factors from the same tables. 

Type 1: u-exchange of an J p = 1/2+ octet baryon 
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Type 3: i-exchangc of a scalar meson (/o) 
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Type 4: i-exchange of a vector meson 
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Type 5: Contact interaction in 7r7V — > pA^ 
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Type 6: t-exchangc of a scalar meson (ao) (scalar coupling) 
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Type 7: N u-exchange in irN — > pN 



(B12) 



V = -ig a — u(p3, A 3 )75^ 9 



+ m-w 



\ 2 - E 1 , - u>2 - ^4 z — Eg — E 1 — E 3 I 
x ^[f (P4,A 4 ) - i ^-^>„ e ;(p 4 ,A 4 )] u (p 1 ,A 1 ) 



(B13) 



Type 8: 7r i-exchange in ttN — > piV 
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Type 9: ^-exchange of a pseudovector meson 
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Type 19: p i-cxchangc in pN —> pN 
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Type 20: A u-exchange in piV — ^ pA^ 
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Type 21: TV u-cxchangc in 7rA — > 7rA 
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Type 22: p i-cxchangc in 7rA — ^ 7rA 
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Type 23: A u-exchange in 7rA — > 7rA 
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Type 25: N u-cxchangc in pN ttA 
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Type 26: ir i-cxchangc in aN —> irA 
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Type 27: N u-cxchangc in aN aN 
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Type 28: a t-exchange in aN —> aN 
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For vertices involving pseudoscalar mesons, we use 



2. Couplings of t- and u-channel exchanges 



olbbp = 0.4 , a P pv = 1 
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The coupling constants of several t- and u-channel ex- 
change diagrams can be related to each other by SU(3) 
flavor symmetry (245| . The coupling of a given meson 
nonet with the octet baryon current 8 b <S> 8b yields two 
products, namely 8s ® 8b <g> 1m and 8b <g> 8 b <B> 8m- The 
coupling to the meson octet 8m depends on two param- 
eters, <7i, <?2, corresponding to the symmetric and anti- 
symmetric representations of 8b <8> 8b- They can be re- 
written in terms of a coupling consta nt q and a mixing 
parameter a in the notation of Rcf. [2451 ] used here as 
well, 
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31 and g 2 can be als o ex pressed in terms of the standard 
D and F couplings j226| : 



D 



30 



40 



■9i 



F 



V6 
24 



92 
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with a = F/{D + F) (definition of a of Ref. [245J). The 
coupling of the baryon octets to the meson singlet 1m 
involves a further independent coupling constant. 



where the index P denotes a pseudoscalar meson and 
V means a vector meson. Wc identify the 77 with the 
octet r/8 and assume the singlet coupling to be zero. This 
leads to a difference of only about 10% in the estimates 
of <7££r/ and g\Ari, assuming an 77 — ?/ mixing angle of 
dp ~ —10°. The corresponding difference in the coupling 
constant gKnK * is o nly about 2%. For the coupling gNNrj, 
following Ref. 246], wc use a phenomenologically small 
value of gNNr] — 0.147 instead of the SU(3) value gNN-q = 



.9b bp 



(4c 



0.343. 



For the coupling of vector mesons to octet baryons, as 
e.g. at the NNp vertex, we use the Lagrangian 



9NNpT 



Jnn p 



2m 



a^d w 



N 



f-p^ (B38) 



which consists of a vector part with 7^ and a tensor part 
with cr A " y . The corresponding coup ling c onst ants are con- 
nected t hroug h Jnnp = gNNp^p [247l |248| ] and we use 



k p = 6.1 H7C 



The couplings of the physical cj and r/> are obtained 
by assuming ideal mixing of the SU(3) singlet- and octet 
states, i.e. of u>\ and cjg. Furthermore, we assume that 
the 4> meson does not couple to the nuclcon (OZI rule) 
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TABLE XI: Couplings for the calculation of the t-/tt-channel 
and contact (cf) diagrams used in this study. The first two 
columns specify the diagram with the exchanged particle (Ex) , 
cf. Figs. [2] and [3] The third column specifies the diagram 
type according to the listing in Sec. IB II the columns 4-6 
show which coupling constants to be used in these expres- 
sions. The values of these couplings are quoted in Sec. IB 21 
and Table IXTTTl The final expressions for the t- and it-channel 
diagrams, given by Eq. (|B3|) . contain also the isospin factors 
(IF), shown here in columns 7 and 8. 



TABLE XII: Continuation of Tabic IXll 
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£ 
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9 
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A 
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n 
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A 
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A 
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A 


14 
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K* 
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Ex 
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TABLE XIII: Couplings not fixed by SU(3) symmetry. The 
values marked with * were obtained in this study by fitting 
to data. For the large size of gNNa in fit B, see the discussion 
in Sec. IIVDI For the sign of g U n P see Ref. [242] . 



Coupling fit A fit B Ref. 



gNNa 

g^Zfa 

g-xpa-1 

1/2 1/2 
g-nr/a gNNa 

1/2 1/2 
g v vfo 9 NNf 

gKKfa 

<?AJVp 

9AAp 



13.85 39.31 * 
9.41 * 
~ gNNp 
10.03 
3.305 
1.336 
16.03 
7.67 



[63,i63j 

w 



[244] 

[122] 
W 



Coupl. fit A fit B Ref. 



gaaa 1.75 5.58 * 

9 aa/ 1.56 1.59 * 

gjVJVai ~ gNN-n [63. 1631 

<W P -10.0 * 

g-nna 1.77 [243] 

<?ppp ~ gwATp [63] 
9aa^ 1.78 [61] 

KAAp 6.1 [61] 



45 



to fix the singlet coupling constant. Note, however, that 
kaon loops prov ide a sm all but non-vanishing effective 
(f>NN coupling [l76l Il77| that is neglected in this ap- 
proach. 

With regard to the actual coupling strengths employed 
in the vector-meson sector, strict SU(3) symmetry is 
broken in various ways. First and foremost, as men- 
tioned before, we do not use genuine p exchange in the 
ttN — > nN channel but determine its contribution from 
correlated tttt exchange. Furthermore, the value of Qnnuj 
as commonly found/needed in studies of the NN inter- 
action exceeds its SU(3) prediction; see the discussions in 
Refs. [M |248| |. Thus, in order to accommodate this as- 
pect, we work with the hypothetical value olbbv — 1-15 
(instead of the standard value ol r hv = 1) based on the ui 
coupling constant given in Ref. |1 TCl| which introduces a 
small SU(3) symmetry breaking in all vector-meson cou- 
plings. Furthermore, we set /nnlu = for simplicity. 

The employed relations for the various coupling con- 
stants are given by the following expressions. Their 
use in the different t- and u-channel expressions of Ap- 
pendix IB II is indicated in the fourth and fifth column of 
Tables [XI] and [XTTJ 

• Couplings for octet baryon, octet baryon and pseu- 



doscalar meson: 


9NN-K 


— gBBP , 


gSNK 


= .9bsp(1 - 2q;_bbp) , 


g\NK 


= -^§-gBBp(l + 2a B BP 




= 2 gBBP a BBP i 




= 7/g gBBP (1 — CKbbp) 




= —gBBP , 


9B.AK 


= 77j gBBp{^at B BP - 1) , 


SAAr, 


= —-J^ gBBp{l — 0>BBP) 


gw v 


= ffSBp(l — O-BBP) , 


with g B BP 


= 0.989 (n3|. 



(B39) 



• Vector coupling for octet baryon, octet baryon and 
vector meson: 

gNNp = gsBv , 

gNNui = 9BBv(4atBBV — 1) ) 
gANK* — — ^jffSBv(l + 2aBw) j 

g^NK* — 5bbv(1 — 2aBBv) , 
ffAAo; = §5-Bsy(5a_Bsy — 2) , 

ffESw = 2gBBVO!BBV , 

= —^gBBvC^O'BBV + 1) j 

= -v / 2.9BBy(2a BB y - 1) , 

flEEp = 2gBBVCtBBV , 

.9AEp = - ttBBv) , (B40) 

with g BB v = 3.25 fl70j 



• Tensor coupling for octet baryon, octet baryon and 
vector meson: 



fNNp = gNNpKp , 

1 a/3 

JaNK* = —l^fNNu) 2~fNNp , 

fsNK* = 2 InNoj + \ fNNp , 

/aAw = |/jVWw — \fNNp , 

/eSlj = \fNNuj + \fNNp , 

/AA0 = — J"/|/aWw ~ -J^fNNp , 

/eE</> = —^fNNu + ^fNNp , 

/eep = \fNNu + \ Jnn p , 

/aE p = —l^pifNN 



2V3 

with Kp = 6.1 |l70l | and Jnnu = 



V3 i 

— jNNp 



(B41) 



• Coupling for pseudoscalar meson, pseudoscalar me- 
son and vector meson: 



g-Kizp 
gxKp 

gKivK" 

gxKui 

gKKrjy 
gKr)K' 

with g PPV = 3.02 [24 



Igppv , 

gppv , 
—gppv , 
gppv , 

V2g PPV , 
-V^gppv 



(B42) 



• Coupling for dccuplct baryon, octet baryon and 
pseudoscalar meson: 



gANn 

gz'NK 

.9E*Ett 

3e*A7t 
gs-AK 

9S*SK 
ffE'Sj) 



with gBBP = 2.13 [17 



— fjDBP j 

= -7S9DBP, 

= -^gosp , 

= -^JDBP , 

= ^5 g DB p , 

= -^goBP, 

= ~^9DBP 



(B43) 



The coupling constants for scalar mesons (/o, <r) are 
not fixed b y SU ( 3) sy mmetry. We take the value of gKK f 
from Refs. |244ll247j . obtained from the hyperon-nucleon 
interaction. The couplings g NN<T , .9<t<t<t, 5ee/ and g A Af„ 
were treated as free parameters and determined in the 
fit. The values can be found in Table IXIIII together 
with other couplings, partly taken from the literature 
and partly fitted in the current study. 
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3. Form factors of t- and u-channel exchanges 



We use 



F(q) = 



A 2 - 



A 2 



(B44) 



for the form factors F of the t- and u-channel exchanges 
in Eq. (|B3[) where m ex is the mass and q the momen- 
tum of the exchanged particle. Powers n = 1,2 provide 
monopole or dipolc form factors. The dipole type applies 
to vertices with p, A, K* , £*, cj and 0. Exceptions are A 
exchange in pN — > pN, where n = 4, and the 7T7T/9 vertex 
in 7rA — > 7rA with p exchange, where n = 1. Otherwise 
the monopole form factor is used. Different powers for 
the form factors are applied to ensure the convergence of 
the integral over the off-shell momenta in the scattering 
equation. 

For the nuclcon w-exchange at the ttNN vertex in 
nN ttN, pN, 7rA, <tN the monopole form factor is 
modified in order to be normalized to unity at the Cheng- 
Dashen point: 



F(q) 



A 2 



' N 



A 2 - ((m 2 N - ml)/m N ) 2 + q 2 
The form factor for the contact terms reads: 



(A 2 - 


m|\ 2 


( A 2 ~ 




\ A 2 4 


-Pi) 




-Pi) 



(B45) 



(B46) 



where P2 (pi) denotes the momentum and m-i {m/Cj the 
mass of the incoming (outgoing) meson. The form factors 
for the correlated ttit exchange potentials can be found 
in Ref. [H. 

Explicit expressions for the form factors of s-channel 
diagrams can be found in Appendix B of Ref. [66| . 

The cut-off values of the form factors for t- and it- 
channel processes, fitted in the present study, are shown 
in Table IXIVl Similar to bare resonance masses and cou- 
plings, we do not attribute any physical meaning to these 
parameters and quote the values only for the sake of re- 
producibility. It is clear that in the present approach 
some values of the cut-offs absorb effects from degrees of 
freedom not explicitly included in the model. If the fit 
prefers a high value for a cut-off, we do not consider the 
corresponding process to be sensitive to high momentum 
physics but rather take it as a sign that some missing 
ingredient of the model is absorbed in that value. 



Appendix C: nN — > nN phase shifts and inelasticities 

In Figs. SS] to US] we show the phase shifts and in- 
elasticities for the reaction ttN — > ttN. The alternative 
representation in terms of partial-wave amplitudes can 
be found in Sec. IIII CI There, one can find the discussion 
on Figs. iH to HI 



S. Dim, Z. Fodor, J. Frison, C. Hoelbling, R. Hoffmann, 
S. D. Katz, S. Krieg, T. Kurth et al., Science 322, 1224 

(2008) . 

C. Alexandrou et al. [ European Twisted Mass Collab- 
oration ], Phys. Rev. D78, 014509 (2008). 

T. Burch, C. Gattringer, L. Y. Glozman, C. Hagen, 

D. Hierl, C. B. Lang and A. Schafer, Phys. Rev. D 74, 
014504 (2006). 

C. Gattringer, C. Hagen, C. B. Lang, M. Limmer, 

D. Mohler and A. Schafer, Phys. Rev. D 79, 054501 

(2009) . 

G. P. Engel, C. B. Lang, M. Limmer, D. Mohler and 

A. Schafer [BGR Collaboration], Phys. Rev. D 82, 
034505 (2010). 

J. Bulava, R. G. Edwards, E. Engelson, B. Joo, H- 
W. Lin, C. Morningstar, D. G. Richards, S. J. Wallace, 
Phys. Rev. D82, 014507 (2010). 

B. J. Menadue, W. Kamleh, D. B. Leinweber, 
M. S. Mahbub, [arXiv: 1109. 67161 [hep-lat]]. 

W. Melnitchouk, S. O. Bilson-Thompson, F. D. R. Bon- 
net, J. N. Hedditch, F. X. Lee, D. B. Leinweber, 
A. G. Williams and J. M. Zanotti et al., Phys. Rev. 
D 67, 114506 (2003). 

N. Mathur, Y. Chen, S. J. Dong, T. Draper, I. Horvath, 
F. X. Lee, K. F. Liu and J. B. Zhang, Phys. Lett. B 
605, 137 (2005). 

M. S. Mahbub et al. [CSSM Lattice Collaboration], 
Phys. Lett. B 707, 389 (2012). 



[11] R. G. Edwards, J. J. Dudek, D. G. Richards and 

S. J. Wallace, Phys. Rev. D 84, 074508 (2011). 
[12] C. D. Roberts and A. G. Williams, Prog. Part. Nucl. 

Phys. 33, 477 (1994). 
[13] C. D. Roberts, Prog. Part. Nucl. Phys. 61, 50 (2008). 
[14] H. L. L. Roberts, L. Chang, I. C. Cloet and 

C. D. Roberts, Few Body Syst. 51, 1 (2011). 
[15] D. J. Wilson, I. C. Cloet, L. Chang and C. D. Roberts, 

Phys. Rev. C 85, 025205 (2012). 
[16] N. Isgur and G. Karl, Phys. Rev. D 18, 4187 (1978). 
[17] S. Capstick and N. Isgur, Phys. Rev. D 34, 2809 (1986). 
[18] U. Loring, B. C. Metsch and H. R. Petry, Eur. Phys. J. 

A 10, 395 (2001). 
[19] M. Ronniger and B. C. Metsch, Eur. Phys. J. A 47, 162 

(2011) . 

[20] F. Gross, G. Ramalho and M. T. Pena, Phys. Rev. D 

85, 093005 (2012). 
[21] G. Ramalho and M. T. Pena, Phys. Rev. D 85, 113014 

(2012) . 

[22] T. Melde, W. Plessas and B. Sengl, Phys. Rev. D 77, 
114002 (2008). 

[23] T. Melde, L. Canton and W. Plessas, Phys. Rev. Lett. 

102, 132002 (2009). 
[24] B. Colli and S. Sirca, Eur. Phys. J. A 47, 61 (2011). 
[25] I. G. Aznauryan and V. D. Burkert, Phys. Rev. C 85, 

055202 (2012). 

[26] N. Kaiser, P. B. Siegel and W. Weise, Phys. Lett. B 
362, 23 (1995). 



47 



TABLE XIV: Cut-offs A for the form factors of the exchange diagrams. The columns Ex specify the exchanged particle and ct 
indicates a contact interaction. 
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